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Introduction 


In Unit 7 you studied a model of vibrating systems which is based upon the concept of 
a perfect spring and leads to the prediction of simple harmonic motion. This unit 
continues the theme of modelling vibration, introducing the three ideas of damping, 
foreing and resonance. 


Damping refers to the inevitable tendency of vibrations to die away whenever a 
system is left to its own devices. This phenomenon was clearly apparent in the 
oscillations of the home-made apparatus of Unit 7 Section 1, but no attempt was made 
there to include damping in the model. In this unit we model damped vibrations by 
including a velocity-dependent resistive force in Newton's second law. 


Forcing is what is needed to keep a system vibrating, in spite of its damping. This can 
arise from a periodic force. For example, a car engine vibrates due to the rotating 
crank even when the car is stationary. A mechanical system can also be forced to 
vibrate by the specific movement of some component of the system. Thus a car on a 
bumpy road is forced to vibrate, at a frequency not of its own choosing, by the bumps 
which its wheels encounter. Similarly, the tone-arm of a record player responds to the 
motion of the stylus in contact with a record. 


Resonance refers to the fact that the amplitude of forced vibrations can be much 
greater at some forcing frequencies than at others. These resonant oscillations can be 
destructive in mechanical systems. Often designers introduce springs and damping 
devices in order to eliminate resonant behaviour over the range of frequencies which a 
piece of machinery is likely to encounter. 


The mathematics required here is that which you studied in Unit 6. In particular, the 
mechanical systems to be considered are modelled by the linear differential equation 


mi +ré+kr = ka, + Pcoswt, 


where m, r, kt», P and w are constants. The case for which r = 0 and P = 0 was 
discussed in Unit 7, and corresponds to simple harmonic motion. Damped vibrations 
are modelled by the addition of the term rz, where r is a positive constant. The 
inclusion of the term P coswt corresponds to forced vibrations. 


Study guide 


The first four sections of this unit should be studied in the order that they appear. The 
main ideas are in Sections 1 and 2. Section 1 uses Newton’s second law, with perfect 
spring forces and resistive forces which are linearly dependent on velocity, to model 
damped vibrations. A periodic applied force is added to this model in Section 2 in order 
to incorporate forced vibrations. This basic model is generalized in two respects in 
Section 3, First, damping forces are discussed which depend on the relative velocity 
between two components of a system. Secondly, we consider forced vibrations which 
arise from the periodic movement of a point of the system other than the particle 
whose motion is sought. 


The audio-tape is associated with Section 3. It will give you practice in using Newton's 
second law to derive the equations of motion which model mechanical systems. The 
television programme is associated with Section 4 and provides a basis for evaluating 
the success of our models of damped and forced vibrations. It will be easier to follow if 
you have already studied the first three sections. 


Finally, Section 5 contains some further exercises on the whole unit. You may use these 
for additional practice, either during the current study week or when you are revising 
for the examination. 
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1 Damped vibrations 


1.1 Introduction 


In Unit 7 we modelled a variety of vibrating systems, and in each case the model led to 
the prediction of simple harmonic motion. But in the real world no system, left to its 
own devices, maintains the oscillations of constant amplitude that are predicted by 
such models; in reality the amplitude decreases with time. In this section we therefore 
introduce a new ingredient to model the observed behaviour of damped vibration. This 
ingredient is a velocity-dependent force. 


One particular system considered in Unit 7 was a home-made apparatus used in a 
simple experiment. The system consisted of a bag of coins suspended from an elastic 
string whose other end was attached to a fixed frame. The bag was pulled down below 
its position of rest and then released, after which it oscillated up and down between 
extreme positions on either side of its equilibrium position. The mathematical model of 
Unit 7, used in an attempt to describe this behaviour, predicts that when the 
displacement of the bag of coins is plotted against time, the resulting curve will be 
sinusoidal, as shown in Figure 1. 


— tA 


a) 


Figure 1 Simple harmonic motion with amplitude A and period >. 


Specifically, the model predicts that the equation of motion is # + w*.x = w*x, and that 
the displacement of the bag is given by + = x, + Acos(wt + @), where x, is the 
equilibrium displacement, A is the amplitude, w is the angular frequency and ¢ is the 
phase angle. 


height of bag + 
above grouni 


static 
height 


- 
time 


Figure 2. The observed behaviour of the home-made apparatus of Unit 7 Section 1, 
showing damped vibration. 
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For this simple harmonic motion both the amplitude and the period r = 27/w are 
constant. How does this prediction tally with observed fact? The experimental results 
shown in Figure 3 of Unit 7 Section 1 supply the answer. The essential features of this 
observed motion are repeated in Figure 2 above. 


In one respect the prediction corresponds with what is observed: the motion does have 
a constant period. But from another point of view the prediction is clearly wrong: the 
motion is not one with a constant amplitude. In fact, the amplitude decreases steadily 
with time and eventually the motion ceases altogether. This observed decrease in 
amplitude is not represented in the simple harmonic model, and in the current section 
we intend to remedy this deficiency. 


1.2 Linear damping 


It is a fair guess that a major reason for the decreasing amplitude of the vibration 
represented in Figure 2 is the effect of air resistance. However, it is known that the 
repeated lengthening and shortening of a rubber band causes an energy loss due to 
internal friction, and this will also contribute to the decrease in amplitude. In the 
context of vibration, frictional effects are commonly referred to as damping and the 
corresponding motion is spoken of as a damped vibration. 


Damping can be due to a number of physical effects and there is a corresponding 
variety of mathematical models. In this unit we concentrate on the resistance 
experienced by a solid body moving through a liquid or gas. As you saw in Unit 4, this 
resistance is commonly modelled as a force opposed to the motion and proportional to 
some power of the speed of the body relative to the fluid. The resistance to the motion 
of the bag of coins through the air is an example of this effect, as is the air resistance 
acting on moving vehicles. 


In Unit 4 two models for air resistance were investigated. They assumed the magnitude 
of the resistive force to be proportional respectively to speed, or to (speed)?. In 
mathematical modelling it is a good principle to start with the simplest available model, 
and not to complicate matters unless this initial model turns out to be inadequate. We 
shall thetefore model the air resistance on the bag of coins as being proportional to the 
speed of the bag (which, since the air is substantially stationary, is also the speed of the 
bag relative to the air). Thus if the velocity of the bag is ¢ then the damping force is 
assumed to have magnitude R = r|z|, where r is a positive constant, and direction 
opposing the motion. This is a much used model called linear damping. It has the 
virtue of simplicity and is mathematically easy to handle, which is such an advantage 
that the model is sometimes used even in situations where it may not be strictly 
appropriate. Any part of the frictional force on the bag of coins which is due to effects 
other than air resistance (such as internal friction in the elastic string) is a case in 
point. In such circumstances one would not expect the model to give results which are 
accuraté in detail, but it is capable of showing the effects of friction in general terms. 


Linear damping 


We use a linear model for damping which stipulates a damping force of 
“magnitude R = r|<|, where « is the velocity of the particle and r is a positive 
constant. The direction of the damping force is opposite to that of the 
velocity. 


1.3. The revised equation of motion 


We now investigate how the inclusion of linear damping alters the predictions derived 
from our model for the oscillating bag of coins, and the first step is to obtain the new 
equation of motion. As in Unit 7, the bag of coins is modelled as a particle of mass m 
and the elastic string as a perfect spring of stiffness k. Here, however, the particle is 
also acted upon by a damping force of magnitude R. Figure 3(a) shows the model 
without specifying any forces. Its purpose is to establish the origin for the displacement 


‘The SI units of the damping 
constant r are Nm~'s (force 
divided by speed), or 
equivalently kgs~'. 
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x of the particle. The most convenient origin for our purposes is the rest position of the 
particle, that is, the position in which it will remain stationary when released. 


mg 
(a) (b) 


Figure 3 


Exercise 1 (revision) 


A particle of mass m is hung from a perfect spring of stiffness k and natural length fo, as 
shown in Figure 3(a). Find the length / of the spring when the particle is in equilibrium. 


{Solution on page 39 | 


As well as choosing the origin at the equilibrium position, the direction of the x-axis is 
taken to be downwards, as indicated in Figure 3(a). This means that: 


(i) x is positive when the particle is below its equilibrium position and negative when 
it is above; 


(ii) % is positive when the particle is moving downwards and negative when it is 
moving upwards; 


(iii) % is positive when the acceleration of the particle is directed downwards and 
negative when the acceleration is directed upwards; 


(iv) the x-component of any force acting on the particle is positive if the force is 
directed downwards and negative if the force is directed upwards. 


Figure 3(b) shows the forces which act on the particle when it is moving downwards 
and the spring is extended. In accordance with the convention used in Units 4 and 7, 
the arrows in Figure 3(b) indicate the directions of the forces and the quantities written 
beside the arrowheads are the force magnitudes. Thus T represents the magnitude of 
the spring force: the spring is in tension, and so it exerts an upward pull on the 
particle. The damping force of magnitude R also acts in the negative .r-direction (since 
its direction is opposite to that of the velocity), and the gravitational force on the 
particle is downwards, as always. Hence, adding the z-components of these forces and 
using Newton's second law, we have 


mg -T — R=mi. (1) 


We next seek to replace the quantities T and R by expressions which involve the 
variables x and é. The result of Exercise 1 shows that when the particle is in 
equilibrium (where x = 0) the spring has extension mg/k. If the particle has position x 
below its equilibrium position then the spring will be further extended by the distance 
x, so that its total extension for this particle position is mg/k + x. The magnitude of 
the spring force will then be the product of its stiffness and extension, that is, 


) =mg+ker. 


Also, since é is positive (with the particle moving downwards), the model of linear 
damping introduced in Subsection 1.2 gives 


R=rlé|=re. 


You tackled this problem as 
Exercise 9 of Unit 7 
Section 2. 
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Substituting these expressions for T and R into Equation (1), we obtain 
mg — (mg + kx) —ré = mi, 
or mé + ra + ka = 0. (2) 


This is the new equation of motion, including the effect of linear damping. Notice that 
all the terms on the left-hand side have the same sign, and that if damping were to be 
removed (by putting r equal to zero) then the remaining equation would be that of the 
simple harmonic model derived in Unit 7. 


Equation (2) was derived by considering a particle moving downwards and an extended 
spring. However, it applies also in other circumstances. Suppose, for example, that the 
spring is still extended but that the particle is moving upwards. Then the situation will 
be as shown in Figure 4 below. 


R 


mg 


Figure 4 


The direction of the spring force in Figure 4 is upwards as before, and the expression 
derived earlier for T in terms of x will be unchanged. The velocity < is now negative, so 
the damping force is directed downwards (in the positive x-direction) and has 
magnitude 
R=r|a|=—rz. 
Thus Equation (1) is replaced by 
mg -~T+R= me, 
but when —ré is substituted for R and mg + kx for T, we have 
mé +ra + kx = 0, 


which is Equation (2) once more. This demonstrates an important feature of the linear 
model for the damping force: the same equation holds throughout the motion, and in 
deriving this equation we need consider only one particular configuration. Usually it is 
simplest to assume that the velocity is positive, or equivalently, that the particle is 
moving in the positive x-direction. Similarly, as with the simple harmonic model in 
Unit 7, there is no loss of generality in assuming that the spring is extended. The case 
when the spring is compressed leads to exactly the same equation of motion. 


Exercise 2 


Redraw Figure 3 for the case in which x measures the displacement of the particle from its 
equilibrium position and the z-axis is directed upwards. For the force diagram, assume that 
the spring is extended and that the particle is moving upwards. Show that the corresponding 
equation of motion for the particle is Equation (2) above. 


[Solution on page 39] 
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The analysis above used the equilibrium position of the particle as the origin of 
position, for the sake of simplicity. It is, of course, possible to use other origins, and it 
is sometimes preferable to measure the displacement of the particle from the fixed end 
of the spring, as shown in Figure 5. You saw a similar alternative for the undamped 
perfect-spring model in Unit 7. 

Example 1 


A particle of mass m is hung from a fixed point by a perfect spring of stiffness k and 
natural length /o. The particle is subjected to a linear damping force with damping 
constant r. 


(i) Find the equation of motion of the particle when its displacement x is measured 
downwards from the fixed end of the spring. 


(ii) Hence find the length /, of the spring when the system is in equilibrium. 


Solution 


(i) The system is illustrated in Figure 5. We assume that the spring is extended and 
the particle is moving downwards, so that the velocity # is positive. (Other cases 
lead to exactly the same equation of motion.) The forces acting on the particle, as 
shown in Figure 6, are: 


(a) the gravitational force, which is directed downwards and has magnitude mg; 


(b) the spring force, which is directed upwards (since the spring is extended) and 
has magnitude T = h(a — Io); 


(c) the damping force, which is directed upwards (since the particle is moving 
downwards) and has magnitude R = ri. 


Newton's second law gives 
mi =mg-—T-—R 
= mg — K(x — lo) — rz, 
or mé +ré+kr = mg + kl, 
which is the required equation of motion. 


(ii) When the particle is in equilibrium, we have ¢ = 0 and # = 0. Substituting these 
into the equation of motion above produces 


ka = mg + klo, 
mg 
a= — +h. 
or r= +10 


Hence the equilibrium length of the spring is |; =lo+mg/k. © 


Example 1(i) shows that when the particle displacement is measured from the fixed end 
of the spring, the corresponding equation of motion is 


me + rt + kx = mg + ko. (3) 


Notice that the left-hand side of this equation is identical in form to that of 
Equation (2). However, Equation (2) is homogeneous whereas Equation (3) has a 
non-zero right-hand side. A particular solution of the inhomogeneous linear differential 
equation (3) is given by 

kay = mg + kilo. 

mg 

or t=, lo, 
which corresponds to the equilibrium position of the particle. This confirms that 
Equations (2) and (3) are alternative descriptions of the same physical behaviour. 


Indeed, if the equation of motion is known with respect to one choice of origin, then we 
can find the equation of motion with respect to any other such choice by a change of 
variable. For example, the equation of motion of our damped vibrating system when 
the displacement of the particle is measured from the fixed end of the spring is given by 


Figure 5 


z 


Figure 6 


mg 
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Equation (3). The displacement of the particle from its equilibrium position, measured 7 
downwards, is (see Figure 7) 
ee ee ing, 
X=2-h=2-(o+7¥). ; I 
From this we have z = X + Ip + mg/k, ¢ = X and # = X. Equation (3) can then be 
m 


written in terms of X as 


mX +rX +k (X+lo+ 72) = mg + ko, 


or mX +rX +kX =0, 


which is just Equation (2) with the variable X in place of x. 


Exercise 3 
‘The equation of motion (Equation (2)) for the damped vibrating system considered above is 
me +ré+ ka =0 4 


when the displacement is measured downwards from the equilibrium position. Let X be the 

upward displacement of the particle from floor level (see Figure 8). Assuming that the 

equilibrium position of the particle is a distance h above the floor, use a change of variable to re 
derive the equation of motion in terms of X. k 


[Solution on page 39 ] 


1.4 The condition for oscillation =nbRmny 


When the bag of coins considered earlier in this section is displaced from its equilibrium 
position and released, it oscillates with a progressive decrease in amplitude, Imagine 
now what would happen if, instead of being suspended in air, the bag were suspended 
in treacle or in thick oil. Intuitively, you can see that its motion would be much more 
sluggish and that vibration might not occur at all. The fact that the system may or WITT? 
may not oscillate, depending on the amount of damping present, is reflected in the 


Figure 8 
types of solution which are obtainable from the equation of motion. i 


The general form of this equation for the model we are considering is 
mé + ra + ka = kare, (4) 


where « is the position of the particle relative to some origin and 2, is the displacement 
of the equilibrium position from that origin. A mechanical system whose equation of 
motion is of this form is called a damped linear harmonic oscillator, or damped 
harmonic oscillator for short. 


In Unit 6 you saw that the general solution of a differential equation such as 
Equation (4) is the sum of any particular solution and the complementary function, 
which is the general solution of the associated homogeneous equation 


me + ra + kr = 0. (5) 


A particular solution of Equation (4) is x = x,, corresponding to the situation in which 
the particle stays at rest in its equilibrium position. To find the condition under which 
the particle may oscillate about this point, it is necessary to solve Equation (5). Indeed 
if x, is taken to be zero, that is, if the origin for x is chosen to be at the equilibrium 
position of the particle, then Equation (5) is itself the equation of motion for the 
system. For simplicity we shall assume for the moment that this is the case, but note 
that the general solution of Equation (4) may be obtained directly from that of 
Equation (5) by addition of the term z-. 


The auxiliary equation corresponding to Equation (5) is Here we are applying 
2 Procedure 1.1 of Unit 6. Note 
mX° +rA+k=0, that this procedure leads to 
an oscillatory solution only in 
case (iii), when the auxiliary 


whose solutions are given by 


=p = ame equation has complex roots. 
rr (6) 
2m 


If the solution for x in terms of time t is to be oscillatory then the roots \ must be 
complex, for which the condition is r? < 4mk. In other words, if the damping constant 
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il 


r is small enough (for example, in air) then the system can vibrate, but if r is too large 


(for example, in treacle) then no oscillation occurs. 


If r? <4mk then Equation (6) can be written as 


ya bivamk =F 
= 2m is 


The complementary function is therefore 


x(t) = sec ( 


= Ae~** cos(Mt + @). 


where 


a > eae 


(7) 


(8) 


and A, @ are constants (with A non-negative and ¢ between —z and 7) whose values 


may be chosen to fit the initial conditions of the system. Equations (7) and (8) give the 


general solution of Equation (5) provided that r? < 4mk. Equation (7) represents an 


oscillatory motion with angular frequency © and a diminishing amplitude Ae~". 
The general shape of its graph is shown in Figure 9, which resembles in all essential 


respects the observed motion shown in Figure 2. Thus the inclusion of damping in the 
model for the motion of a bag of coins has enabled us to predict an important feature 
of the motion which did not show up in the simple harmonic model, The revised model 


is therefore an improvement on the first one. 


Figure 9 


The period of the damped motion is tr = 27/9, with © given by Equations (8). The 
corresponding period for undamped motion was found in Unit 7 Section 2 to be 27/wo, 


where wo = \/k/m is the angular frequency of simple harmonic motion. Not 


surprisingly, the expression (8) for 2 reduces to \/k/m when the damping constant r is 


equated to zero. 


Exercise 4 


Show that, for a given mass m and spring stiffness k, the period of linearly damped oscillations 


is greater than the period of undamped oscillations. 
[Solution on page 39] 


The amplitude of the oscillations described by Equation (7) decreases steadily due to 
the exponential factor e~?*, and the greater the damping (that is, the larger the value 


of r), the faster this amplitude decreases. To verify this, we shall look at the ratio 


a(t +7) 
a(t) * 


According to Procedure 1,1 
of Unit 6, the right-hand side 
of Equation (7) would have 
been written as 

e"(Beos Mt + Csin Mt), 
where p, 9 are as in 
Equations}(8) and B, C are 
arbitrary constants. 
However, we have written the 
term in brackets here in the 
alternative form given by 
Equation (5) of Unit 7 
Section 2. 


In the remainder of this unit 
we use wo to denote the 
angular frequency of simple 
harmonic motion. The 
symbol w, which was 
employed for this purpose in 
Unit 7, will be given a 
different meaning in 

Section 2. 
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where 7 = 27/9 is the period of damped oscillations. Using Equation (7), we have 
a(t +7) = Ae?" cos (Q(t + 7) + @) 
= Ae") cos( Nt + 2x + 9). 
Now cos(Qt + 27 + @) = cos(Nt + @), and so 


attr) _ Aevnttr) 
a(t) Ae~#t 
As p > 0 and r > 0, we have 


SieTetttr=t) = errr, (9) 


ef <1 and 2(t+7) < x(t). 


Equation (9) gives e~?” as the constant ratio of any two non-zero values of «x one period 
apart; it applies, for example, to two consecutive maximum displacements on the same 
side of the rest position. As the damping constant r increases, it can be seen from 
Equations (8) that p increases while 2 decreases. Hence 7 increases and e~?* decreases, 
confirming that increased damping produces a faster decrease in amplitude. The box 
below summarizes the results of this subsection. 


Linearly damped vibrations 

The equation of motion of a damped harmonic oscillator, namely 
mé +ré + kar = 0, 

has an oscillating solution provided that 
r? < Amk. 

The general solution in this case is 
x(t) = Ae~*' cos(Mt + 0), 


where 


2m 


and A, @ are arbitrary constants. All solutions of this form may be obtained 
with A non-negative and @ between —7 and 7. 


The angular frequency of these damped vibrations is 2, and their period is 
T = 27/Q. The ratio of two non-zero values of x one period apart is 
a(t+7) _ 


=ePr 
a(t) , 


We next return for a final look at the experimental results obtained from the suspended 
bag of coins, showing how measurements of the amplitude can be used to estimate a 
value for the damping constant r. 
With 20 coins in the bag, it was found that the maximum displacement of the bag from 
its equilibrium position decreased by a factor of about 10 over 8 complete cycles. 
Taking x to represent displacement from the equilibrium position, this means that 
a(t +87) 
San = 0.1, 
Now the left-hand side of this equation can be re-expressed as 
a(t+8r)  a(t+7r) a(t+r) 
x Kix = 
a(t+7r) x(t+6r) a(t) 


Thus we have e~*?7 ~ 0.1, from which 


(e7°7)® = ef, 


pr = —L log, 0.1 = 0.288. 
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Using Equations (8) and r = 27/9 gives 

Lar 20 x 2m 
Eee 

a 0.2887 

d4mk—r? ~ 4x? * 

This equation may be solved to obtain a value for r provided that numerical values for 

the constants m and k are available. In Unit 7 Subsection 2.4 it was stated that each 

coin has a mass of 0.01 kg, so 20 coins have a mass of 0.2kg. It was shown also that the 

elastic string can be modelled by a perfect spring of stiffness 2.28N m~'. Solution of 

the equation above with m = 0.2 and k = 2.28 results in the estimate 


r ~0.062Nm~'s 


= 0.288, 


or 


for the damping constant. 


As you saw above, the condition for oscillation is r? < 4mk. For the current case, this 
means that oscillation is possible provided that 


r?<4x02%2.28, or r<1.35. 


The estimate 0.062 for r is less than 3; of this upper bound, showing that quite a small 
damping force is enough to make the amplitude decrease fairly rapidly. The effect of 
small damping on the angular frequency and hence on the period is less marked, as you 
are asked to demonstrate in the next exercise. 


Exercise 5 

Estimate the angular frequency of oscillations when there are 20 coins in the bag, as predicted 
by 

(i) the model for linearly damped motion, using the value of r derived above; 

(ii) the model for undamped motion developed in Unit 7. 


Exercise 6 


If as before m = 0.2, k = 2.28 and r = 0.062, find the values for A and @ in Equation (7) which 
correspond to the initial conditions (0) = —0.05 and #(0) =0. (Recall that A should be 
non-negative and that ¢ lies between —x and =.) 


Exercise 7 


Figure 10 shows a highly simplified model of a car. The particle of mass m represents the car 
body. This is supported by a perfect spring of stiffness k and a ‘damper’, which together 
represent the suspension. The contribution of the wheels is neglected and only the vertical car body 
motion is modelled. The damper contains a piston which moves through oil inside the 

cylinder. This causes a resistance to the motion of the piston which is approximately 

proportional to its velocity relative to the cylinder, the constant of proportionality being r. 


As part of a quick test of the suspension by a potential buyer when the car is stationary on a damper} 


horizontal road, the body is depressed 0.1 m below its rest position and then released. If 

m = 800, k = 5 x 10* and r =9 10°, all in the appropriate SI units, establish whether the 
ensuing motion is oscillatory and, if so, find the distance of the car body below its rest position Figure 10 
after one complete cycle. 


{Solutions on page 39 | 


1.5 Weak, critical and strong damping 


In Subsection 1.3 we derived the equation of motion 
ma +ré+ ka =0 


for a mechanical system consisting of a particle of mass m moving under the influence 
of a perfect spring of stiffness k and linear damping with damping constant r, where 
the displacement of the particle is measured from its equilibrium position. In 
Subsection 1.4 the general solution of this equation was found for the case in which the 
inequality r? < 4mk is satisfied. There are, of course, solutions also for the cases when 
r? =4mk and r? > 4mk, and you are asked to derive these in the next exercise. 
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Exercise 8 

Find the general solution of the differential equation 
me +ra + kr =0, 

where m, r and k are positive constants, when 

(i) 1? <4mk, 

(ii) 7? =4mk, 

(iii) r? > 4mk. 


Solution on page 40] 


The results of this exercise are summarized below. 


Weak, critical and strong damping 
The equation of motion 
mé+ré+kx =0 
has the following solution. 
(i) If? <4mk (weak damping) then the general solution is 
x(t) = Ae™** cos(Mt + ¢). 


where 


els 
e* mm 
Here A and ¢ are constants; all solutions of this form may be obtained 
with A non-negative and ¢ between —7 and 7. 
(ii) If? = 4mk (critical damping) then the general solution is 
x(t) = Be" + Cte", 
where p = r/2m and B, C are arbitrary constants. 
(iii) If r? > 4mk (strong damping) then the general solution is 
a(t) = Be?#* + Ce7P2, 
where B, C are arbitrary constants and 


b r— Vr? —4mk mi r+ Vr? —4mk 
= A = , 
2m 2m 


The case r? < 4mk, which is called weak damping, was discussed in the previous 
subsection. The general solution 
x(t) = Ae" cos(Qt + ¢) 


is the product of an exponential decay and a sinusoidal oscillation of angular frequency 
Q, So the particle oscillates backwards and forwards about its equilibrium position with 
period 7 = 27/2, much as it would when there is no damping, but with the amplitude 
of the oscillations decaying exponentially. The graphs of typical undamped and 
weakly-damped oscillations are shown in Figure 11 below. 


When r? > 4mk, which is called strong damping, the solution is 
a(t) = Be-Pt 4 Ce-Pat, 


where p; and po are positive constants. This solution is the sum of two decaying 
exponential terms. So from any initial conditions the particle returns to its equilibrium 
position without oscillation; the damping force dominates, and prevents any vibration. 
A graph for the strongly-damped case appears in Figure 11 below. 


In the case r? = 4mk, which is called critical damping, the solution is 
x(t) = Be" + Cte", 


Part (i) of this exercise is 
revision of Subsection 1,4. 
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where p is a positive constant. Both parts of this solution represent a decaying motion 
and so the motion again fails to be oscillatory, as shown in Figure 11. Given any initial 
displacement and velocity, the particle returns asymptotically to its equilibrium 
position just as for strong damping. As suggested by the figure, critical damping gives 
the fastest return to the equilibrium position. Hence many physical devices, such as 
recoil mechanisms for large artillery guns, mechanical suspensions, door-closing 
mechanisms and weighing scales, are designed to have critical or near-critical damping. 


critically 


at 


Figure 11 


Exercise 9 


A gum has a barrel, of mass 10° kg, which starts to recoil at a speed of 30ms™' when the gun 
is held horizontally and fired. Immediately after firing, the barrel starts to compress a spring 
which brings the barrel instantaneously to rest after it has recoiled a distance of 2m. In this 
phase of the barrel’s motion no damping is present. A linear damping mechanism is then 
engaged to ensure that the barrel returns to its initial position without overshooting the mark, 
(i) What is the spring stiffness? 

(ii) What damping constant is required for critical damping? 

{Solution on page 40 | 


Summary of Section 1 


1, Damping can be modelled by a damping force of magnitude r|i|, where r is a 
positive constant and # is the velocity of the particle, and direction opposed to the 
velocity of the particle. The SI units of the damping constant are Nm~'s (or 
equivalently kgs~'). 


2. The equation of motion of a particle of mass m moving under the influence of a 
perfect spring of stiffness / and linear damping with damping constant r is 


me +ré + kx =0, 


where «x is the displacement of the particle from its equilibrium position. Such a 
system is called a damped harmonic oscillator. 


3. For weak damping (r? < 4mk) the general solution of the equation of motion is 
a(t) = Ae? cos (Nt + 4). 


where 


_ vin 

2m ‘ 

All solutions of this form may be obtained with the arbitrary constants A and @ 
restricted so that A is non-negative and @ lies between —z and 7. This motion is 
oscillatory with period t = 27/9. The amplitude of the oscillations decreases 
exponentially, so that 


a(t +7)/x(t) = e7*7. 
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4. For critical damping (r? = 4mk) the general solution is 
a(t) = Be" + Cte, 


where p = r/(2m) and B, C are arbitrary constants. This case gives the fastest 
return to the equilibrium position. 


5. For strong damping (r* > 4mk) the general solution is 
x(t) = Be7?* + Ce7P2", 


where B, C are arbitrary constants and 


r— vr? —4mk r+vr? —4mk 
nS 2m bp eS 2m : 


2 Forced vibrations 


2.1 Introduction 


So far in this unit we have looked at the behaviour of a damped harmonic oscillator 
which moves freely after being started off from an initial position and with an initial 
velocity. We now go on to consider the motion of such a system when it is maintained 
in a state of vibration by some external physical agency. This ‘agency’ is modelled by a 
force changing continuously with time. The model of a car considered in Exercise 7 of 
Section 1 suggests one example: if the engine is started with the car stationary then it 
exerts an alternating force on its supports and hence on the car body, which vibrates 
for as long as the engine is running. As another example, think of a voltmeter which is 
used to measure a voltage varying with time. As the voltage increases and decreases, 
the pointer of the instrument will swing backwards and forwards. In both of these 
cases, the external force applied to the vibrating mass may vary with time in quite a 
complicated manner. We shall assume that the forcing agent is periodic and in that 
case, as will be shown later in this course (Unit 31), the force can be expressed as the 
sum of a number of sinusoidal terms. It follows that if we know the effect of a 
sinusoidal force we can, at least in principle, cope with any periodic variation by 
superimposing the effects of as many sinusoidal components as necessary, In this 
section, therefore, we shall concentrate on the motion of the simple damped model 
when it is subjected to a sinusoidal force. 


2.2 The equation of motion 


Post 
mg 
(a) (b) 


Figure 1 


This ‘superposition of effects’ 
is dealt with by applying the 
superposition principle of 
Unit 6 Subsection 2.4 to 
solve the differential equation 
involved. 
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Figure (a) shows the model system of Figure 3 of Section 1, where the origin for 
position x is chosen to be at the particle’s equilibrium position and the direction of the 
w-axis is taken to be downwards. This system is as before acted upon by the force due 
to gravity, a perfect spring and linear damping, but now there is in addition an applied 
sinusoidal force with x-component Pcoswt. The positive constant P is the amplitude 
(maximum magnitude) of the applied force, and w is its angular frequency. These four 
forces are shown in Figure 1(b) for the case in which the spring is extended, the particle 
is moving downwards and P coswt is positive. However, the same equation of motion 
results also for the other possible configurations of the system. 


As before, the magnitudes of the spring and damping forces for the case shown are 
respectively T = mg + kx and R =r. Application of Newton’s second law to Figure 1 
therefore gives the equation of motion 


mg + Pcoswt — rz — (mg + kx) = mi, 
or mi +ré+ kr = Pcosut. 


qa) 


As you may recall from Unit 6, the general solution of this equation is the sum of two 
parts: the complementary function, which is the general solution of the associated 
homogeneous equation (with a zero right-hand side), and a particular solution of the 
original equation, which depends upon the form of this equation’s right-hand side. The 
complementary function for Equation (1) was found in Section 1; in this section we 
concentrate on the particular solution. 

When, as here, the equation models a physical system (with m, r and k positive), the 
amplitude of the complementary function decreases steadily and eventually becomes 
negligibly small. For this reason the complementary function is often called the 
transient part of the solution of Equation (1). On the other hand, the particular 
solution is persistent since, as you will see, it has a constant amplitude. It therefore 
remains in evidence long after the transient part of the solution has become negligible, 
which explains why the particular solution is known as the steady-state part of the 
solution. 


Since the right-hand side of Equation (1) is sinusoidal, the easiest way to find its 
steady-state solution is to use the phasor method, This is based upon Euler’s formula, 


e'” = cos0 + isind. 
Putting wt in place of @ gives 
elt = coswt + isinwt, 


so that coswt is the real part of e*. It follows that Pcoswt = Re(Pe™*), A particular 
solution of Equation (1) is now obtained by trying a solution of the form 
2 = Beoswt + C sinwt, or equivalently, 

wx = Acos(wt + ¢), 
This is sinusoidal, of the same frequency as the external force but out of phase with it 
by an angle @. Application of Euler's formula with wt + ¢ in place of 0 leads to 

x = Re(Ae'**)) — Re( Ae'*e™*)_ 


We next put Ae’? = Z, where Z is a complex constant (the phasor of the sinusoidal 
solution). The modulus of Z is A, the amplitude of the steady-state solution. The 
argument of Z is @, the phase of the steady-state solution. For x and its first two 
derivatives we then have 

x = Re(Ze“*), 

i = Re(iwZe*), 

# = Re(i*w*Ze™*) = Re(—w*Ze™"). 


(2) 


Substitution of these expressions into Equation (1) produces 
Re{(—mw® + irw + k)Ze"*] = Re( Pe™*). 


Note our notation for the 
three types of angular 
frequency occurring in this 
unit: (i) wo for simple 
harmonic motion (undamped 
and unforced); (ii) 2 for 
damped harmonic motion 
(unforced); (iii) w for an 
applied force (and for the 
steady-state oscillation which 
it causes). 


See the box on page 14, 


This behaviour was 
illustfated in Figure 11 of 
Section 1. 


See Procedure 2.2(a) in 
Unit 6. 


Unit 5 Subsection 4.1 
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Hence Equation (2) will represent a solution of Equation (1) if 
(—mw? + irw + k)Z = P, 


giving 

= Pe oa P 

~ =mw? +ire+k (k— mu?) + ire 

_ Pl(k — mw?) — irw} 
ae a (k = mw?)? + 1r2w? * 
This is the phasor of the steady-state solution for x. Its relationship to P (the phasor 
of the applied force) is shown graphically in Figure 2. It may be expressed as Figure 2 
Z=u-iv, where 

ent 
ae P(k — mw") ae Prw 


~ (k= mu?) ra? 


Hence the modulus and argument of Z are respectively 


Gamat ra 


A=|2|= V+ = (3) 
and @ = Arg(Z) = —arccos (a) = —arecos (4) 
The steady-state solution of Equation (1) is therefore 

2 = Acos(wt +4), (5) 


where A and @ are given by Equations (3) and (4) respectively. This sinusoidal 
solution, with amplitude A, is the ‘output’ from the system corresponding to the 
‘input’ supplied by the forcing term with amplitude P. The angular frequency w of the 
output is the same as that of the input, but between input and output there is a time 
delay of —/w, as illustrated in Figure 3 below. 


me | 


Figure 3 


Exercise 1 

Derive Equations (3)-(5) without the use of phasors, by substituting the trial solution 
«= Beoswt + Csinwt 

into Equation (1), (‘This is the alternative method of Procedure 2.2 in Unit 6.) 
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Exercise 2 


Consider once more the model car described in Exercise 7 of Section 1, for which m = 800, 

k =5 x 10* and r = 9 x 10°, all in the appropriate SI units. While the car remains stationary 
on a horizontal road, its running engine produces a force on the car body which may be 
considered as sinusoidal, with amplitude 800N and frequency 2.5 Hz, (Recall from Unit 7 
Subsection 2.3 that 1 Hz = 1 cycle per second.) 


(i) What is the corresponding value of the angular frequency w in rads~!? 
What is the amplitude of the steady-state vibration experienced by the car body? 


(iii) What is the time lag between the maximum value of the engine force being exerted in the 
upwards direction and the car body next reaching its maximum displacement in that 
direction? 

Exercise 3 

The equation of motion of a mechanical system is 

2% + Se + 14x = cost. 
Use the phasor method to find the steady-state oscillations of this system in the form 
a(t) = Acos(wt + ¢). 


[Solutions on page 40 | 


You have seen above how the phasor method can be used to find the steady-state 
vibrations of a model system when the external force has 2-component P coswt. The 
method employs the relationship 


P coswt = Re( Pe"). 


A similar phasor approach can be adopted when the applied force is of the more general 
form P cos(wt + €), or of the equivalent form M coswt + N sinwt. In these cases we use 
the equations 
Pcos(wt +) = Re(Pe'*e'*"), 
Mcoswt + Nsinwt = Re ((M —iN)e’**) 
respectively. You are asked to consider the second of these possibilities in the exercise 
below, 
Exercise 4 
(i) Verify that 
Re ((M —iN)e™') = Mcoswt + Nsinwt. 
(ii) Use the phasor method to find the steady-state solution of the equation of motion 
#+ 2% + 5a = 2cos2t + 9sin 2t, 
giving your answer in the form 
x(t) = Mcoswt + Nsinwt. 
What are the amplitude and phase of this vibration? 
[Solution on page 41] 


2.3. The damping ratio 


In the previous subsection we used the phasor method to show that the equation of 
motion of the forced and damped harmonic oscillator, 


mit +ré+ke = Pcoswt, (1) 
has the steady-state solution 
a = Acos(wt + 4), (5) 
where 
P 
2————— 6 
Vik ma?) + ra? (6) 


k— mu? 
(R= mut rar 


and ¢@ = —arccos ( (7) 
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The applied force here is completely described by two parameters, its amplitude P and 
angular frequency w. For a particular set of values of m, r and k, the steady-state 
solution %r is a function of both P and w, and we wish now to examine the nature of 
this functional dependence. This amounts to investigating how the output parameters 
A and ¢ depend upon the input parameters P and w. As has already been noted, the 
angular frequency w of the output is identical to that of the input. 


Notice first that the phase angle @ is independent of the input amplitude P while, for a 
fixed value of the input frequency w, the output amplitude A is proportional to P. It is 
not so straightforward to see the effect upon A and @ of varying either the angular 
frequency w or the amount of damping present, but this is the principal aim of the 
current subsection. 


We proceed by reorganizing the information contained in Equations (6) and (7), so that 
the quantities m, r and k are replaced by combinations which permit a clearer view of 
how A and ¢ react to changes in w. In order to provide some motivation for the choice 
of these ‘combinations’ of m, r and k, we return first to Equation (1). Dividing this 
equation by m gives 


#+7e+2= F cosust. (8) 


The coefficient of x is now k/m, a quantity encountered previously. In Unit 7 
Subsection 2.3 you saw that the angular frequency of the unforced and undamped 
perfect-spring system was wo = 80 we may replace k/m by w2. In this context 
we shall refer to wy as the undamped angular frequency. 


The next step is to rewrite the damping term. You saw in Section 1 that the 
complementary function of the equation of motion is oscillatory if and only if the 
damping constant r is less than an upper limit given by r? = 4mk, or r = 2Vmk. This 
leads us to define the damping ratio a of the system by 
r 

Oink” 
The damping constant r provides an absolute value for the damping force per unit 
speed exerted on the particle. On the other hand, the damping ratio a gives a measure 
of the importance of damping relative to the mass m and spring stiffness k of the 
system. Note that a is a dimensionless quantity, since the dimensions of r (Nm~'s, or 
kgs~!) are identical to those of mk. 


From the way in which a is defined, and from the results of Subsection 1.5, it follows 
that a system moving in accordance with the complementary function of Equation (8) is 


(i) weakly damped (oscillatory) if a < 1, 
(ii) critically damped if a = 1, 
(iii) strongly damped if « > 1. 


Figure 4 shows graphs of the complementary function for a particular system for four 
different values of a; the initial conditions are the same for all of them. This is 
essentially a repetition of Figure 11 in Section 1. 


Th 


a 


Figure 4 


As mentioned before, we use 
wo in this unit to denote the 
natural (unforced, 
undamped) angular 
frequency, rather than w as 
in Unit 7, since currently w 
represents the angular 
frequency of the forcing term. 
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From the equations 


r 
= 1 id = , 9 
so m ce & 2Vmk Ka 
we have 
r _ 2aVvmk 


m m 
so that Equation (8) can now be written as 


» er ee 2 
# + Lawyt + wor = — cosut. (10) 


Exercise 5 

Determine the phasor for the steady-state solution of Equation (10). Hence find expressions 
for the amplitude and phase angle of this motion. 

Exercise 6 


A linearly damped perfect-spring system with sinusoidal applied force has the equation of 
motion 


Ai + 9% + 100x = 9 cos 6. 
(i) Find the values of the undamped angular frequency and the damping ratio. 
(ii) Determine the amplitude, phase and angular frequency of the steady-state solution, 
[Solutions on page 41} 


In Exercise 5 you showed that the steady-state solution of Equation (10) is 
a = Acos(wt + @), where 


Pim 


—— ———— uy 
(wa — w?)? + daPuzu? aut) 
ow} —w? 
and =a SS | 12 
an = —arecos ( Kat — atte ER) (12) 


These expressions could also have been obtained directly from Equations (6) and (7), 
using Equations (9) to replace k by mw? and r by 2amwo. 


The forced and damped harmonic oscillator 
The equation of motion 
mi +ra+ kx = Pcoswt 


of the forced and damped harmonic oscillator can be rewritten as 
- Ree 
#+ 2awod + wor = fa coswt, 


in terms of the undamped angular frequency 


k 
wo = 4/ — 
m 


and the damping ratio 
r 


a= 


2Vmk° 
The steady-state solution x = Acos(wt + @) of this system has amplitude 
a) 
Vei— ow) + 4atagat 
and phase 


w2 —w? 


aa ( Vu? — 0)? + 40%a2u? 
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The best way to show the significance of Equations (11) and (12) is to draw graphs. In 
the case of Equation (12), what we want to demonstrate is the effect on the value of @ 
of varying w and a. But the equation also contains wo, so one possibility would be to 
plot, for each of a range of values of wo, a set of curves of @ against w, each separate 
curve corresponding to a particular value of a. This might mean drawing a large 
number of sets of curves, which would require much calculation and plotting. In order 
to simplify the procedure, we note that Equation (12) can be rewritten as 


~ cm ( 


which shows that @ can be regarded as a function of the two variables a and w/w. We 
therefore plot graphs of @ against w/w for particular values of a, as shown in 

Figure 5(b) below. This enables all of the required information to be presented by 
drawing just one set of curves. 


The same sort of approach can be applied to Equation (11), except that the task of 
representing graphically the variation of A with w and a is further complicated here by 
the presence of m and P as well as wy. Noting that 

PEE 2P 


=>— =47, 
m km oF 


and putting M = Ak/P, Equation (11) can be rewritten as 


a (13) 


V= w/o? ? + 407(0/is0)? 
where, as before, the right-hand side is a function of the two variables a and w/w. The 
resulting graphs of M against w/wo for various values of a are shown in Figure 5(a) 
below. 


Notice that M, like a and w/wo, is a dimensionless quantity. For a given spring stiffness 
k, the ratio of the output displacement amplitude A to the input force amplitude P can 
be found from the value of M. 


Ms 


(a) 


(b) 


Figure 5 
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Exercise 7 


Calculate the value of M which corresponds to the steady-state solution of the differential 
equation given in Exercise 6 above. 


[Solution on page 42] 


2.4 Resonance 


The most noticeable feature of the graph of M against w/wo in Figure 5(a) is that for 
low values of a it shows a definite peak near w/wo = 1. This means that, when the 
angular frequency w of the sinusoidal external force is close to the undamped angular 
frequency wo of the spring-mass system, the steady-state amplitude of the displacement 
of the particle becomes large compared with its values for the same force amplitude 
over most of the frequency range. The motion of the system at these relatively large 
amplitudes is known as resonance. As a increases, the resonance peak becomes less 
pronounced and moves towards lower values of w/w. When a = 1/2, the peak is at 
w/wy = 0 (as you are asked to show in Exercise 10 below). For a > 1/2 the curve has 
negative slope for all positive values of w/w. 


Resonance is found in many physical systems and may be considered desirable, as in 
the case of a radio receiver which is tuned to an input signal of a particular carrier 
frequency, so that inputs at other frequencies are neglected. However, there are other 
situations where resonance may be most undesirable. In the case of a car, large vertical 
movements of the car body will certainly be uncomfortable and may also be dangerous 
and destructive, 


Figure 5(a) may be used to represent the response of an instrument which measures 
and records a physical quantity such as, for example, electrical voltage. The movement 
of the instrument's pointer (and attached parts of the mechanism) can be modelled by 
the movement of a particle, and the voltage causes, and is proportional to, the external 
applied force. If the voltages to be measured are periodic and have a range of different 
frequencies, then it is preferable to avoid re-calibrating the instrument for every 
different frequency; indeed, it might be the case that the voltage consists of several 
sinusoidal functions of different frequencies superimposed. Hence the same value of M 
is required for as wide a frequency range as possible. Of the curves drawn in 

Figure 5(a), that for @ = 0.6 meets most closely the condition that M should be 
approximately constant over the greatest possible frequency range. The value usually 
chosen in practice is @ = 0.65. In order to achieve this damping ratio, the instrument 
movements often incorporate deliberately added linear, or nearly linear, damping of a 
mechanical or electrical nature. 


In order for the output to be a faithful representation of the input we need to have the 
same value of M for all the frequency components, as already explained, so that each 
component of the input makes the appropriate contribution to the output. In addition, 
we need to ensure that the time shift between different components (the interval 
between the instants at which different components reach their maxima, say) is not 
affected by passing through the instrument. If it were, the output would be a distorted 
version of the input. So, to avoid distortion, one wants the time lag introduced by the 
instrument to be the same for each component. This time lag equals —¢/w (see 

Figure 3, page 18), so that we want the phase ¢ to be directly proportional to the input 
frequency w. As it happens, for a = 0.65 this is almost true over a frequency range 
similar to that for which M is approximately constant. Figure 5(b) illustrates this. 


Figure 5(b) also shows that, for all values of a and w/wo, the phase ¢ lies between —7 
and 0, and that @ = —7/2 when w = wy. 


Exercise 8 

A voltmeter is modelled by a forced and damped harmonic oscillator with damping ratio 

@ = 0.6 and undamped angular frequency wo = 20rads~*. Estimate from Figure 5(a) the 
frequency range over which the instrumert will reproduce the amplitude of a sinusoidal input 
voltage with a relative error of less than 10%. 


Exercise 9 


Use the graphs in Figure 5(a) to check your answers to parts (ii) and (iii) of Exercise 2 
(page 19). 
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Exercise 10 


Write the right-hand side of Equation (13) in terms of a and p = w/wo. Using this expression 
show that, for a given value of a, the maximum value of M occurs 


(i) when p= V1 — 202, if a < 1/V2; 
(ii) when p=0, if a > 1/2. 
[Solutions on page 42} 


Summary of Section 2 
1. The equation of motion of a damped harmonic oscillator subject to an applied 
sinusoidal force of amplitude P and angular frequency w can be written as 
me + ki +rx = Pcoswt. 


The steady-state forced vibrations of this system can be found by using the 
phasor method. 


2. This equation of motion can be rewritten in the form 
. BE 
# + 2awpk + wpx = — coswt, 
m 


where the undamped angular frequency o is defined (as in Unit 7) by 


Je 
wo = —. 
m 


and the damping ratio a is defined by 
r 


2Vimk- 


3. The steady-state vibration of the forced and damped harmonic oscillator is given 
by 


a= 


a« = Acos(wt + ¢), 


where 
ae P/m 
Vw) = w?)?? + 402g? 
wi — w? 
and = = —arccos | ————* ——______ } . 
Viet) data? 


4, For a <1/¥2, the amplitude of the forced vibrations exhibit a maximum at a 
particular value of the forcing angular frequency w. This effect is known as 
resonance. 
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3 The perfect dashpot 


3.1 Introducing the dashpot 


In Section 1 we augmented the Unit 7 model of a particle-spring system by the 
inclusion of linear damping, and in Section 2 an external force acting upon the particle 
was added. This section shows how the linear model for damping can be applied to 
configurations of particle, spring and damping which differ from those encountered 
previously. These configurations arise from situations in which an external force acts 
upon a point of the system other than the particle itself. 


For example, consider the vertical oscillations experienced by a motor car as a result of 
its motion along an undulating road, This may be modelled by considering the car to 
be at rest with the road surface moving beneath it, the effect of the undulations being 
represented by a sinusoidal displacement of the wheels about the mean road level. In 
predicting the vibrations of the car body under these circumstances we need to include 
a linear damping force which is proportional to the relative velocity between the body 
and wheels, rather than a force proportional to the absolute vertical velocity as in the 
earlier model for a stationary car. Also, we now have to incorporate an external force 
described by the vertical displacement of the wheels rather than a force acting directly 
upon the car body, 


The main purpose of this section is to enable you to write down the correct equation of 
motion for the particle when the system consists of any configuration involving a 
perfect spring, linear damping and an applied force. As a first step, Figure 1 shows a 
diagrammatic means of indicating linear damping. It represents a piston inside a 
cylinder, or dashpot. One form of dashpot is used as the shock absorber on a car (see 
Exercise 7 in Section 1), where the cylinder is full of oil and the relative motion of the 
piston forces the oil to flow through the annular space between piston and cylinder. 
This produces a resistance to the relative motion. For modelling purposes we assume 
that damping is produced by a perfect dashpot, for which the resistance is 
proportional to the relative velocity of the dashpot’s ends (or equivalently, to the rate of 
change of the dashpot’s length). The constant of proportionality is denoted as before 
by r, but is now referred to as the dashpot constant rather than the damping 
constant. If the instantaneous length of the perfect dashpot is / then the magnitude of 
the damping force produced is r//\. 


The occurrence of a perfect dashpot on a diagram indicates the presence of linear 
damping, irrespective of whether the damping being modelled is produced by an actual 
dashpot or by some other physical means. In this respect the role of the perfect 
dashpot is analogous to that of the perfect spring, which indicates the presence of a 
force depending linearly on the relative displacement of two end-points, whether or not 
the real force being modelled is applied by the agency of a spring. 


Figure 2 shows how we could draw a perfect dashpot alongside a perfect spring to 
represent the model of Section 1 for the suspended bag of coins. The damping effects 
here are due not to a physical dashpot but to air resistance and internal friction in the 
rubber string. Since the point of suspension is fixed, the damping force magnitude r|i| 
in this instance is equal to r|#|, where x is the displacement of the particle from any 
fixed point. However, this will not always be the case. 


Figure 3 shows a perfect dashpot with both ends in motion. The positions of these ends 
with respect to some fixed point are denoted by x and y, so that their velocities are 
respectively ¢ and yj, Since the dashpot length is | = x — y, the magnitude of the 
damping force exerted by the dashpot is r|i| = r|i — g|. Now if ¢ > g, that is, if the 
dashpot is lengthening, then the force exerted on whatever is connected to the dashpot 
at either end is directed towards the dashpot. (Imagine that you are pulling the two 
ends of the dashpot apart with your hands. The dashpot will resist by exerting forces 
trying to stop your hands moving apart. The force on your right hand is to the left, and 
the force on your left hand is to the right.) By the same token, if ¢ < j, that is, if the 
dashpot is shortening, then the forces on objects attached to the ends will be directed 
away from the dashpot, again resisting any change to the length of the dashpot. 


Figure 1 
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Example 1 
a 
oO ra 
——e m 
i 
1 r r 
r— f 
y ' 
r 
Figure 4 


A particle is attached to one end of a perfect dashpot, as shown in Figure 4. The 
displacements of the particle and of the other end of the dashpot from a fixed point O 
are denoted respectively by x and y, where x > y. Find the x-component of the 
dashpot force acting on the particle when 


(i) the length of the dashpot is increasing! 


(ii) the length of the dashpot is decreasing. 


Solution 
The length of the dashpot is 1 = x —y, and sol =#—y. 


(i) When the length of the dashpot is increasing, / is positive and so |I| =/1 =z — 9. 
The magnitude of the dashpot force is therefore R = r|i| = r(¢ — i). In this case 
the dashpot force is directed towards the dashpot, as shown in Figure 5, trying to 
resist the increase in the length of the dashpot. Hence the 2-component of the 
dashpot force is 


—R=-r(#— 4). 


(ii) When the length of the dashpot is decreasing, / is negative and so 
\i| = -i = —(# — y). Hence the magnitude of the dashpot force is 
R= ri = —r(é—). Again, this force tries to resist the change in the length of 
the dashpot, and so is directed away from the dashpot, as indicated in Figure 6. 
‘Therefore the x-component of the dashpot force is 


R=-r(é-y). O 


This example demonstrates that the -component of the dashpot force takes the same 
form whether the length of the dashpot is increasing or decreasing. Consequently, the 
equation of motion of a particle attached to a dashpot will also have the same form in 
these two cases. We shall therefore usually assume, when deriving an equation of 
motion, that the dashpot length is increasing (so that the dashpot force is directed 
towards the centre of the dashpot). However, the resulting equation of motion will be 
valid also for the case when the length of the dashpot is decreasing. 


The perfect dashpot 


The perfect dashpot is a diagrammatic representation of a type of force, 
rather than of the mechanism which produces it. It represents a resistive force 
proportional to the rate of change of the dashpot’s length /. In other words, 
the force has magnitude R = r|l|, where r is a positive constant known as the 
dashpot constant. 


If the dashpot’s length is increasing (so that / is positive), the dashpot force is 
directed towards the centre of the dashpot. If the length of the dashpot is 
decreasing (so that / is negative), the force is directed away from the centre of 
the dashpot. In either case the dashpot force resists any change in the length 
of the dashpot. 


The equation of motion for a particle attached to a perfect dashpot can be 
derived by assuming that the dashpot is increasing in length. 


m 


Figure 6 
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Example 2 


Figure 7 


A particle of mass m is attached to a fixed point O by a horizontal perfect spring of 
natural length [o and stiffness k. The particle is also attached by a horizontal dashpot 
of dashpot constant r to a point whose displacement with respect to the fixed point O 
is y, as shown in Figure 7. Find the equation of motion for the particle, assuming that 
it is otherwise free to move in a horizontal direction. 


Solution 


In deriving the equation of motion we shall assume that the spring is extended and that 
the length of the dashpot is increasing, so that the forces act in the directions indicated 
in Figure 7. We shall measure the displacement xr of the particle from the fixed point 
O, The extension of the spring is x — lo, so the magnitude of the spring force is 


T =k(x —lo). 


‘The length of the dashpot is I = x — y, and so / = + — j. Now [ is positive, since the 
length of the dashpot is increasing. Hence the magnitude of the dashpot force is 


R=ri=r(é-¥). 


By Newton's second law, the equation of motion for the particle is 


mé=-T-R 

= —h(x — Io) —r(z — 9). 
or = mé+ré+khe=htry, O 
Exercise 1 


Find the equation of motion for the particle in Figure 2, taking x to be the downward 
displacement of the particle from the fixed ends of the spring and dashpot, and Ip to be the 
natural length of the spring. Compare your answer with Equation (3) in Section 1 (page 9). 


Exercise 2 
ype ft 
0 1 = 
® { 
Pooswt 
k 
Figure 8 


Find the equation of motion for the particle in Figure 8, which is attached to a perfect dashpot 
and perfect spring (of natural length lo) as well as being subjected to an external force whose 
r-component is P coswt. 


[Solutions on page 42] 


Consider once more the equation of motion derived in Example 2 for the situation 
shown in Figure 7. Note that, as it stands, this equation of motion cannot be solved to 
obtain an expression for the function x(t), since the ry term on the right-hand side is an 
unknown function of t. This situation can be altered by specifying the position function 
y(t) of the left-hand end of the dashpot. This left-hand end is referred to as the forcing 
point of the system, since it is the point at which an external force is applied. 
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There are essentially three different configurations for a single-particle system, 
connected to one perfect spring and one perfect dashpot, which has an external force 
applied at a point other than the particle: 


(i) the external force may be applied to the other end of the spring, with the other 
end of the dashpot held fixed; 


(ii) the external force may be applied to the other end of the dashpot, with the other 
end of the spring held fixed; 


(iii) the other ends of the spring and dashpot may be joined together, and the external 
force applied to this joint end. 


Example 2 looked at an instance of case (ii). We next consider an occurrence of 

case (iii), including specification of the forcing point motion. Provided that this motion 
is sinusoidal, the mathematics required to solve the equation of motion for the particle 
is much the same as that employed earlier. 


Figure 9 is a diagram of the highly simplified model of a car which was mentioned at 
the beginning of this section. This model is obtained by neglecting the mass of the 
wheels and concentrating the wheel suspensions into one pair of perfect spring and 
perfect dashpot elements. The spring and the dashpot are assumed to have the same 
length at all times, with their lower ends being at the same horizontal level. The ‘input’ 
y(t) is a sinusoidal displacement of the common lower end of the spring and dashpot. 
This input represents the undulations of the road, which the car is traversing at a 
certain speed, The ‘output’ x(t) is the vertical displacement of the particle of mass m 
which represents the car body. Both x and y are taken to be positive upwards and 
measured from a fixed horizontal level, which is the mean height of the road surface. 


Assuming that the spring is extended and that the dashpot is increasing in length, the 
directions of the forces on the car body are as shown in Figure 9. The extension of the 
spring is 2 — y — Ip, and the rate of increase of the dashpot’s length is # — j. So the 
magnitudes of the spring and dashpot forces are respectively 


T=k(x-y-lo) and R=r(it—y). 
The gravitational force, of magnitude mg, acts vertically downwards. The equation of 
motion of the car body is therefore 
mg—-T—R 
mg — K(x —y—lh) —r(é - 4), 


me = 


or me + ra + ka = klo — mg + ky + ry. () 
If the undulations of the road are sinusoidal about its mean level and the car travels at 
constant speed, then the forcing point motion will also be sinusoidal. If these forcing 
point oscillations have amplitude Ap and angular frequency w then y = Ap coswt, giving 

mé + ra + ka = klo — mg + Ao(kcoswt — rw sinwt). (2) 
This equation for forced and damped motion is similar to those which you met in 
Section 2. Its steady-state solution has the form 

a(t) = Ip — mg/k + Acos(wt + ¢), 
where Ig — ng/k is the height at which the particle could remain in equilibrium if no 
external force were applied. This constant part of the solution is obtained by finding a 
particular solution of the equation 

mé +ré + ka = kly — ma. 
The sinusoidal part of the solution is derived by finding a particular solution of the 
equation 

mé +raé + ka = Ao(kcoswt — rwsinwt). (3) 


Exercise 3 


Derive an expression for the phasor of the steady-state solution of Equation (3), and hence 
derive expressions for the amplitude A and phase angle ¢ of the forced vibrations described by 
Equation (2). 


[Solution on page 42] 


O| W777, 


Figure 9 
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Under certain circumstances the position x of the particle relative to some fixed point 
may be of less interest than its position = relative to the forcing point. Since > = x — y, 
it is a simple matter to obtain the appropriate differential equation for = directly from 
that for x. For example, the height of the car body above the road surface in the model 
of Figure 9 is = = x —y. On replacing x by = + y in Equation (1), we obtain 


mi +ré+kz = kly — mg — mj. 


Given that y = Ag coswt, this equation for =(t) can be solved in the same manner as 
was adopted for Equation (2). 


3.2 Deriving equations of motion (Audio-tape Subsection) 
The audio-tape provides further practice in setting up equations of motion, analysing 
two model spring-dashpot systems which are slightly more complicated than those 


considered previously. 


Start the audio-tape when you are ready. 


Note that during the tape the 
symbol v (Greek ‘nu’) is used 
in place of w for the angular 

frequency of the forcing term, 


30 MST204 8.3 


@ First System 
particle 


one of mass m 
spring - sinusoidal 


stiffness k eS external 
force 


1 


two dashpots - Horizontal motion only 
constants r, & tp es 


Equation of motion ? ~- What forces act ? <- Variables ? 


6 System at Rest 


Take x as particle position. 
rest position of particle 
natural ; origin for x 


length of direction 
spring for x: 


nn 


What are 
the forces 
© magnitude ? 
© direction? 


(4 ) The Force Diagram 
0! x 


x 
ne 
ae nx ne force 

lashpo 
ok 


Equation of motion ? | 
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(5) The Equation of Motion 


x-components of forces are Newton's 2nd law: 


@ external: Pcos vt mx = x-component 
© spring: -kx of total force 
© dashpots: ty Xt X 


| mx = Pcos vt =ke =r X-1.x | 


Re-arrange : mx +(r,+r,) x +kx = Pcos vt 


or | mx + rx +kx =P cos vt (r=ry +r) 


hens | rest position of 


particle: origin for x 


Forces balance: 
mg=T=k (L)- to) 
S>L= fo +mg/k. 


Frame 5 equation applies again 


8) Second Second System —at Rest at Rest 


two springs, Springs are at 
one dashpot their natural lengths 


9] In Motion 


Assume x > 0, 
x>0, y>x 


What are 
the force 


© spring | © spring 2 © dashpot 
extension = fe) extension = eal length change rate = (| 
tension = (=r tension = | force magnitude = 
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® The Force Magnitudes 


extension=| X extension 
© spring 1 a a © spring 2 Ea 


fension = fension = 


Hache a length change rate= |_x Diagram of 
© dashpo' : 
P force magnitude = | rx forces ? 


® The Force Diagram 


: 


kx<eeH 


Spring forces 
Equation of motion ? 


® The Equation of Motion 


x-components of forces are 


Newton's 2nd law: 
mx = x—component 
© spring |: ~kyx of total force 
© spring 2: +kz (y-x) 


mx = ~1ri-k, x +k2 (y-x) 


Re-arrange: mx + rx +(k,+k,)x=k2y 


or | mk +rxt+kx=kzy (k=k, +k, ) 


© dashpot: -rx 
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The remaining exercise in this section again asks you to derive an equation of motion, 
and illustrates how the response graphs of Figure 5 in Section 2 (page 22) can be 
adapted for use in different circumstances. 


Exercise 4 


Figure 10 shows a model of a solid body being moved vertically through water. The perfect 4 
dashpot represents the resistance of the water, and the perfect spring represents the resilience 

of a towing rope, (Note that a perfect spring can reasonably be used as a model for a rope 

only if the rope is in tension, since a rope, unlike a spring, cannot sustain a compressive force.) 


(i) Initially the top end of the spring is held at a fixed position. If the natural length of the r 
spring is lo, what is its length when the particle is in equilibrium? k 


(ii) The top end of the spring is now moved along a vertical straight line, with displacement y 
measured upwards from its previous fixed position. The corresponding displacement x of 
the particle is measured upwards from its equilibrium position, By making suitable 
assumptions, draw the diagram of forces acting on the particle, and derive the equation 
of motion for the particle. 
si Figure 10 
(iii) Suppose that k = 200Nm~' and m = 4kg, and that the damping ratio is a = 0.2. 
Suppose further that the top end of the spring is given a sinusoidal displacement about 
its original fixed position, of amplitude 0.05 m and frequency 2Hz. By comparing the 
equation of motion derived in part (ii) with Equation (1) of Section 2 (page 17), use 
Figure 5(a) in Section 2 (page 22) to estimate the amplitude of the particle’s steady-state 
oscillations. 


{Solution on page 43 | 


Summary of Section 3 


1. A perfect dashpot represents a resistive force which is proportional to the 
relative velocity between two components of a mechanical system. 


2. The dashpot force has magnitude R = r|/|, where | is the length of the dashpot 
and r is a positive constant known as the dashpot constant. If the dashpot’s 
length is increasing, then the dashpot force is directed towards the centre of the 
dashpot; if the length is decreasing, then the force is directed away from the centre 
of the dashpot. 


3. The equation of motion of a particle attached to a perfect dashpot can be derived 
by assuming that the dashpot’s length is increasing. This equation of motion 
applies also to the case where the dashpot is contracting. 


4 Off the record: resonance and 
damping (Television Section) 


Now watch the television programme. 

Read the following after viewing the programme. CC): 
8 

Allan Solomon introduced the programme by saying that to understand why the 

tone-arm of a record player can misbehave, we need to understand damping and forcing TVv8 


in simpler mechanical systems. Thereafter the programme divided naturally into three 
parts, each of which considered the vibrations of a particular apparatus: 


Part 1: the unforced, damped vibrations of an apparatus modelled by Figure 1; 
Part 2: the forced vibrations of an apparatus modelled by Figure 2; z 
Part 3: the forced vertical vibrations of a tone-arm, modelled by Figure 3. 


Part 1 of the programme showed how a simple mass-spring-dashpot system “ 
behaved, and introduced the model for it shown in Figure 1. David Broadhurst found m 
that a dashpot containing air made little difference to the oscillations observed without iT 


a dashpot: the motion still resembled simple harmonic motion (see Figure 4(a)). With 

water in the dashpot the damping was more noticeable, though the motion was still 

oscillatory (Figure 4(b)). Strong damping occurred with oil in the dashpot, when there 

were no oscillations but just a simple decay to the equilibrium position (Figure 4(c)). Figure 1 
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(a) air 


(b) water 


eins 


{c) oil 


Figure 4 Output traces from a simple mass-spring-dashpot system, 


It was then shown that the equation of motion for the model in Figure 1 was 
mi tra + kx = mg +klo, 


where Io is the natural length of the spring. This equation was solved, and graphs of 
the predicted motion were shown for three values of the damping ratio a = r/(2Vmk) 
(see Figure 5, in which to = lo + mg/k is the equilibrium position of the particle). For 
a <1 damped oscillations were forecast, while for a > 1 an absence of oscillations was 
predicted. Critical damping, with a = 1, gave the fastest decay to the static position. 
The observed results were clearly in qualitative agreement with the model's predictions. 
One quantitative test is the subject of the exercise below. 


Exercise 1 


In Subsection 1.4 it was shown that, for the model of Figure 1, the ratio of particle 
displacements from the equilibrium position at times one period 7 apart is 


a(t+7) ~er 
2° 
where p = r/(2m). It follows from this that 


a(t+nr) _ enter 
a(t) 
Thus if Ao, Ai, A2,-.. are the heights above the mean of successive peaks on the graph of x(t), 


and Aj, Aj, A3,... are the depths below the mean of successive troughs (with Aj between Ao 
and Aj, say), then 


An = Ape"? and Al, = Age™"*7, 


for n = 1,2. 


from which 
An + Aj, = (Ao + Ape”, 
Ao+ Ao \ _ 
= a BER 


Using the data provided below for the water-filled dashpot, test this prediction of the model 
by drawing an appropriate graph. 


n SESEAES Ee 
100(An + Ay) | 5.1 | 2.2 | 1.3 | 0.8 | 0.6 


(Solution on page 44 | 


m 


Figure 3 


This is Equation (3) of 
Section 1, which you derived 
once more in Exercise 1 of 
Section 3. 


ro 


O 12 8 4 5 6 Two 


Figure 5 

Predicted motions for the 
model system shown in 
Figure 1. 


‘We add the two previous 
equations here because it is 
easier to measure Ay + / 
from the trace than it is to 
measure A, or Aj, 


individually. 
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Part 2 of the programme dealt with the forced oscillations of an apparatus 
modelled by Figure 2. The design of this apparatus is indicated in Figure 6. 


pens i ‘pinged arms 


output 


input 


Figure 6 The apparatus used in Part 2 of the television programme. 


The motor-driven cam made the bottom of the spring move up and down sinusoidally 
with a fixed amplitude and at a forcing frequency that could be varied. This was the 
input oscillation, 


Y = Yo + Ao coswt, 
recorded by the lower of the two traces. The steady-state output oscillation, 
t= x9 + Acos(wt +), 


was recorded by the upper trace. In addition to the spring connecting the two arms of 
the apparatus, there was a dashpot connecting the upper arm with a point at a fixed 
level. The geometry is such that the effective damping constant can be increased by 
sliding the dashpot to the left. 


Initially the output was studied for three different forcing frequencies, with the same 
amount of damping. The first frequency gave an amplitude ratio A/Ao greater than 
unity. At a higher frequency the amplitude ratio was even greater, but at a higher 
frequency still it decreased. Thus the system exhibited resonance, since the amplitude 
ratio had a maximum greater than unity (see Figure 7). 


T2Ay 


input VLD PP DAP ADD PDP PDP DP PP DLP PIMA 


Figure 7 Input and output traces for the apparatus of Figure 6 with light damping, 
near the natural frequency. 


The apparatus was modelled by the system of Figure 2, leading to the equation of 
motion 
mi + ra +ka = klo —mg + ky, (1) 
where /p is the natural length of the spring. With y = yo + Ap coswt, Equation (1) 
becomes 
mi +ré + kx = k(Io + yo — mg/k) + kAg cos wt. (2) 
This has steady-state solution 
© = 29 + Acos(wt + 9), 
where ro = lo + yo — mg/k and 
A i 1 
= -- =: (3) 
Ao (I= (w/uo)?)? + (2aw/wo) Figure 8 
Graphs of the amplitude ratio from Equation (3) were provided for values of the sae predictions of 


tion (3) for different 
damping ratio a between 0.1 and 1.0 (see Figure 8). These showed that with Eee) oie 


sufficiently large damping, an absence of resonance was predicted. ratio a. 
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Exercise 2 
(i) Show that the equation of motion for the system in Figure 2 is given by Equation (1). 


(ii) By comparison with the analysis carried out in Section 2, show that the sinusoidal term 
in the steady-state solution of Equation (2) has the amplitude ratio described by 
Equation (3). 

(iii) State a lower bound for the values of the damping ratio a which correspond to 
non-resonant behaviour in this model system. 


[Solution on page 44 | 


Data taken from the apparatus of Figure 6 are shown in Figure 9 below. It was found 
that a reasonable fit to these data could be obtained with the graph for a = 0.2 in 
Figure 8. The prediction of no resonance for sufficiently large damping was also tested. 
With the dashpot moved considerably to the left, a sweep of input frequencies revealed 
a steady decrease of output amplitude with increasing frequency. This indicated an 
absence of resonance, in qualitative agreement with the prediction. 


A/Ag } 
a4 


0 10 2 » 0 w 


Figure 9 Amplitude ratio data obtained from the apparatus of Figure 6. 


Part 3 of the programme dealt with the forced vertical oscillations of the tone-arm 
of a record player. It was shown that the large demonstration apparatus of Figure 10 
exhibited damped vibrations when unforced and resonant vibrations when forced. Like 
the previous apparatus, this one involved the pivoting motion of an arm. 


hinged arm 


springy connection 


motor-driven 
7 forcing point 


Figure 10 Demonstration apparatus for Part 3 of the programme. 
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The motion of the tone-arm was modelled by the system shown in Figure 3, in which 
the connection between the stylus and tone-arm is assumed to involve a perfect dashpot 
(linear damping) as well as a perfect spring. In the following exercises you are asked to 
find the equation of motion of this system, and the amplitude of the steady-state 
output when it is subjected to a sinusoidal input vibration. Note that interest centres 
on the motion of the tone-arm relative to the record surface, as described by the 
variable = = x — y. 


Exercise 3 
(i) Show that the equation of motion of the model system shown in Figure 3 is 
me + ré + kx = klo — mg +ry + ky. 
(ii) Hence show that the displacement = = x — y of the particle relative to the forcing point 
satisfies the differential equation . 
mi +ré +kz = kl — mg — mj. 


Exercise 4 


If the system shown in Figure 3 is subjected to a sinusoidal input vibration y = Ao coswt, 
show that the amplitude ratio for steady-state vibrations of the particle relative to the forcing 
point is given by 

A (w/wo)* 


Ao = (w/a)??? + 40" wn) 


where wo = \/k/m is the undamped angular frequency and a = r/(2V’mk) is the damping 
ratio, 


[Solutions on page 44} 


In the previous exercise you derived an expression which predicts the amplitude ratio 
for the tone-arm. Some corresponding graphs are shown in Figure 11. It was stated 
that the vertical oscillations of a real tone-arm can be reasonably accounted for by this 
model with a value of about 0.2 for a. This is true, at low frequencies, but there are a 
couple of remarks which need to be made. First, the cartridge detects two sorts of 
motion: relative vertical motion between the stylus and tone-arm, and relative 
horizontal motion (the existence of these two signals is what makes stereo reproduction 
possible). Secondly, the arm is not perfectly rigid and can therefore flex in response to 
the motion of the stylus (this is important at frequencies higher than those we consider). 
‘These two factors mean that the real system has more ‘degrees of freedom’ than have 
been incorporated by modelling the arm as a particle moving in one dimension. 


Despite these limitations, the model has certain satisfactory features. It can be shown 
that it predicts resonance for a < 1/\/2 (using the method of Exercise 10 in Section 2, 
page 24), and it suggests that we should increase the damping if we want to reduce the 
amplitude ratio of the system near the natural frequency. However, this natural 
frequency is typically 60 radians per second (that is, 60/(27) ~ 10 oscillations per 
second, or 10 hertz). The frequencies at which vibrations are forced by the undulations 
in the record grooves, producing the music that we hear, are between ten and a 
thousand times greater in value. This led to the question of what could be causing the 
record-player to misbehave so badly. 


The fault turned out to be a bad warp on the record. All records are warped to some 
extent, causing potentially an unwanted resonant response to forcing at low frequencies. 
This effect can be reduced by increasing the damping. It was shown that a similar 
tone-arm on a second turntable did not exhibit such large forced vibrations. The 
difference between the tone-arms was that the second one had a specially designed 
miniature dashpot connecting the arm to the surface of the record. (To make the 
comparison a fair one, a sliding weight on the second arm had been adjusted to 
compensate for the additional weight of the dashpot.) 


oO 1 2 3 4u/uy 


Figure 11 
The predictions from the 
result of Exercise 4, 
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5 End of unit exercises 


Exercise 1 


A small heavy object is hung from the ceiling by a spring, and it is found that the spring is 
extended by 19.6cm in the static position. The object is then attached to the floor by two 
identical dashpots, chosen to give a total effective damping ratio a = 1. 


(i) Devise a model for this system, and obtain its equation of motion, giving numerical 
values for the coefficients. Take the value of g to be 9.8ms~?. 


(ii) One of the two dashpots is now removed. What is the new equation of motion? Predict 
the period of the damped vibrations that can now be observed. If the object is released 
from rest with the spring at its natural length, what is the maximum subsequent spring 
extension predicted by the model? 

Exercise 2 


Suppose that, in addition to the forces acting on the object in Exercise 1(ii), there is an 
externally-produced magnetic force, acting downwards with magnitude 


Po + Pcoswt (Po > P > 0). 
Show that the resulting forced vibrations have amplitude 
d 
=P ye +P 


where p = w/wo and d is the distance by which a constant downward force of magnitude P 
extends the spring. 


Exercise 3 


Figure 1 shows a model of a seismometer, which consists of a particle of mass m suspended 
inside a rigid box by a perfect spring and attached also to a perfect dashpot. The box is given 
a vertical sinusoidal motion of angular frequency w and amplitude Apo, and the displacement = 
of the particle relative to the box is recorded. This displacement is zero when the particle is in 
equilibrium. 


(i) Obtain the equation of motion of the particle relative to the box. 
(ii) Derive the steady-state solution of this equation of motion. 


(iii) Show that, if the damping ratio a and the ratio wo/w of undamped to forcing angular 
frequencies are both small, then = ~ —y. (This result demonstrates that the device can 
be used to record the effect of an earthquake, provided that the damping is light and the 
undamped angular frequency is small compared with the frequencies of ground 
movement.) 


Exercise 4 
This exercise again refers to the model of Figure 1, and to the results of Exercise 3. 


(i) Show that, if the damping ratio a and the ratio w/wo of forcing to undamped angular 
frequencies are both small, then the relative displacement = of the particle is 
approximately proportional to the acceleration of the box, with the constant of 
proportionality being independent of the forcing frequency. (This result demonstrates 
that the device can be used as an accelerometer, provided that the damping is light and 
the undamped angular frequency is large compared with the frequencies of accelerated 
motion to be measured.) 


(ii) Ifa =w/wo = 0.1, calculate the ratio of the approximate amplitude of the box 
acceleration (as indicated by the instrument) to the true amplitude. 


[Solutions on page 45] 
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Appendix: Solutions to the exercises 


Solutions to the exercises in Section 1 


js 


Figure 1 


In the equilibrium position, the spring will be extended and 
so the spring force will be upwards. The magnitude of this 
force is T = k(l; — lo). The only other force on the particle is 
that due to gravity, As the particle is in equilibrium we have 
T=mg, 
which leads to 
4, = lo +:mg/k. 
2. The figure required is shown below. The force diagram 


assumes, as stated in the question, that the spring is 
extended and the particle is moving upwards. 


mg 
Figure 2 


As in Exercise 1 we have T = mg when the particle is at its 
equilibrium position. If the particle is below this position 
then, since x < 0, the additional tension in the spring is 
k(—w) =—kx. Alternatively, if the particle is above its 
equilibrium position (but with the spring still extended) 
then, since x > 0, the decrease in tension is kx. In either 
case, the total tension in the spring is given by 

T = mg — ke. 
‘The particle is moving upwards, so that z > 0 and 

R=r\i|=re. 
Newton's second law then gives 

mi =T -—R—mg 

= ka —ré 


or mi+ré+kx=0, 
which is Equation (2) of Section 1 once more. 


3. The particle is at a distance x below the equilibrium 

position, which is at a height h above the floor. So the height 

of the particle above the floor is X = h —-r. Hence we have 
z=h-X, *=-X 2 

Substituting these expressions into the equation 
mi + rt + kr =0, 


and £=-. 


we obtaii 
m(—X)+7r(—X) +h(h-X)=0 
or mX +rX +kX =kh. 


4. The angular frequency 2 of damped oscillations is 


= Vimk=? _ 
= 


where wy = \/k/m is the undamped angular frequency. 
Hence 2? < wé, so that r = 27/2 is greater than the 
undamped period 27/wo. 


5. (i) The damped angular frequency is 


Vamk —r? 
n= — —.. 
2m 
Here m = 0.2, k = 2.28 and r = 0.062, so that 


(il) The undamped angular frequency is 


=e a en - 
wo = yf = 1/ TS = 3.38radse. 


6. Equation (7) of Section 1 is 
a(t) = Ae~* cos(Mt +). 
‘The derivative of this is 
2(t) = —pAe™" cos(Mt + 6) — NAe~"* sin(Nt + d). 
Using the given initial conditions, we have 
2(0) = Acos¢ = —0.05, (e) 
#(0) = —pAcos@ — QAsing = 0. (2) 
Now Q ~ 3.37 from Exercise 5(i) and, using the first of 
Equations (8) in Section 1, 


r _ 0.062 
P= an = page = Mls. 


Use of Equation (1) to substitute for Acos@ in Equation (2) 
therefore leads to 
ing — 2:05 _ 0.05 x 0.155 _ 

Asing = =>? ~ A = 0.0023. (3) 
Equations (1) and (3) can be solved by the method of Unit 7 
Subsection 2.3, yielding 

A= /(-0.05)? + (0.0023)? ~ 0.050, 

@ = arccos(—0.05/A) ~ 3.1. 


7. The motion will be oscillatory if r? < 4mk. Here we have 


r=9x10°, m=800, k=5~x 104, 
so that 
r?=8.1x 10", 


Amk =4 x 800 x 5 x 10" = 1.6 x 10°. 


Hence the inequality is satisfied, and the motion is 
oscillatory. 
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The ratio of displacements one period apart is 
a(t+T) _ pr 


x(t) 
where 
r Qn Vamk — 1? 
Rg ee 
Therefore 


ee eee 2n x9 x 10* 
Or Vimk= JT6x 10° — 8.1% 10" 


From this it follows that 
Cmagire Bee ae 


Hence one complete cycle after being released the 
displacement of the car body below its rest position will be 


0.1 x 1.73 x 107% = 1.73 x 10-* m, 
or about 0.2mm. 


8. For the differential equation 
mi + rt + kr = 0, 
the auxiliary equation is 
mM +rA\+hk=0 
with roots 
ya te Rank 
2m 


In each of the following three cases we apply Procedure 1.1 
of Unit 6. 


(i) r? <4mk (complex roots for A) 
The roots are \ = ~p + i, where 
7? 


tin 
ai Se) an 


so the general solution of the differential equation is 
a(t) = Ae~" cos(Mt + %). 
Here A and @ are arbitrary constants, though all solutions of 


this form may be obtained with A non-negative and @ 
between —zx and 7. This case is called weak damping. 


(ii)? = 4mk (equal roots for A) 


The auxiliary equation has just one solution, A = —p, where 
p=r/(2m). So the differential equation has general solution 


a(t) = Be" + Cte", 


where B and C are arbitrary constants. This case is called 
critical damping. 


(iii) r? > 4mk (distinct real roots for \) 
The roots are \ = —p; and \ = —p2, where 
vr? —4mk r+ Vr? —4mk 


2m 2m 
Note that both p; and p2 are positive. The general solution is 


a(t) = Be + Ce", 
where B and C are arbitrary constants. This case is called 
strong damping. 


a= and p2= 


9. (i) The most efficient way of finding the stiffness k is 
to use energy conservation (Unit 7 Section 3), since there is 
no damping in operation between the firing of the gun and 
maximum spring compression. Equating the initial kinetic 
energy, with m = 10° and v = 30, to the potential energy for 
4 compression of 2 metres, we obtain 
aie ana 

$x 10° x 30 = $k x 2? 

or k=2.25 x 10°Nm™. 


ii) For critical damping we require r? = 4mk, that is, 
J req 


r= 2Vmk = 2y/10° x 2.25 x 105 = 3 x 10°. 


‘Thus a damping constant of 3 x 10*Nm™'s is required. 
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Solutions to the exercises in Section 2 


1. We look for a particular (steady-state) solution of the 
form 


a= Beosut + Csinwt, 
whose first two derivatives are 
#=—Businwt + Cwcoswt, 
= —Bu* coswt — Cw" sinut. 
Substitution of these expressions into the differential 
equation 
mi +ré+kx = Peoswt 
produces 
m(—Bw” coswt — Cw* sinwt) 
+7r(Cwcoswt — Bu sinwt) 
+ k(Booswt + Csinwt) = Pcoswt. 


We can now compare the coefficients of the cosine and sine 
terms in turn on both sides of this equation: the coefficient 
of coswt on the left-hand side must equal P on the right, and 
the coefficient of sinwt on the left-hand side must equal zero 
(since there is no sinwt term on the right). This gives 


—mBu? +rCu +kB = P, a) 

—mCu? — rBw+kC =0. (2) 
From Equation (2) we have (k —mw*)C = rwB, or 

,_ _TwB 

(ad (3) 


Substituting for C in Equation (1), we obtain 
E oe 
(k — mw’ )B+ 7B =P 
which leads to 
pe ee mw?) 
(k = mu?) + 2? * 
It follows from Equation (3) that 
Prw 
w= (k= mw?)? + rw?” 
Having found the steady-state solution in the form 
x= Beoswt + Csinwt, 
it remains to translate this into the form 
t= Acos(wt +6). 
where, from Unit 7 Subsection 2.3, 
A= /B+C?, 
and = -—arceos(B/A) (since C > 0). 
Using the expressions above for 3 and C, we find 
P 


V (k= mw?)? + 12? 


and a 


A= 


mu? 
Vk = ma??? +20? } 
in agreement with Equations (3) and (4) of Section 1. 


{Notice that this approach takes longer than the phasor 
method.| 


2. (i) The angular frequency is 
xf (where f is frequency in Hz) 
Qn x 2.5 = 15.7rads', 


(ii) From Equation (3) of Section 2, the amplitude is 
A 


With m = 800, k = 5 x 104, r=9 x 10%, P = 800 and 
w ~ 15.7, we have 


k— mw? ~ -1.47 x 10° 
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and 
o_O 
(LAT x 10°)? + (9 x 10° x 15.7)? 
Fe et pees 3 
= Vaarx ioe ~ Z0ax toe = 992% 10m 


so the amplitude of the vibration is about 4mm. 


(iii) From Equation (4) of Section 2, the phase angle of the 
vibration is 


( k— mu? ) 
= —arccos | —————- —— ] . 
Vk = ma??? + ru? 
Using the values above, we obtain 

=1L4T x 10° 
2.04 x 10° 
—arccos(—0.722) 
—2.38rad (or —136"). 
‘The time lag is therefore 
b 2.38 
wo 157 


6 = anes ( 


015s. 


3. ‘The applied force is cost = Re(e“'), with angular 
frequency w = 1. We look for a steady-state solution of the 
form wr = Re(Ze""), where Z is a complex constant (namely, 
the phasor of xr). Successive differentiations of this give 


ut =Re(iZe") and #=Re(—Ze"). 
Substitution of these expressions into the equation of motion 
produces ~ 


Re(—2Ze'' + 5iZe"' + 14Ze"") = Re(e') 
or Re((12 +5i)Ze"*) = Re(e"). 
‘This equation will be satisfied provided that 


(12-4 5i)Z =1, 
for which we require 
as Et eee er 
12+ Ie Te ses 


— arceos( #2 
-0.395 rad (or ~22.6°). 
Hence the steady-state solution is 
+ =Re(te'*e") =Re(de''**) 
= jcos(t + 4), 
where @ ~ —0.395 rad. 


"we have 


4. (i) Using Euler's formula to expand e 
(M =iN)e“* 
=(M ~iN)(coswt + isinwt) 
=(M coswt + Nsinwt) + i(M sinwt — N coswt), 
so that 
Re((M —iN)e“*) = M coswt + N sinwt. 
(ii) The applied force is (using the above) 
2cos 2t + 9sin 2t = Re((2 — 9i)e?"*). 
‘The angular frequency of this applied force is w = 2, so we 
look for a steady-state solution of the form x = Re(Ze”"), 
whose first two derivatives are 


#=Re(2iZe™") and # = Re(—4Ze*"). 
Substituting these expressions into the equation of motion, 
we have 


Re(—4Ze™"* + 4iZe™* + 5Ze™"*) = Re((2 — 9i)e™") 
or Re((1+4i)Ze™"*) = Re((2 — 9i)e™*). 
This equation will be satisfied if 

(14 4i)Z = (2-91), 
for which we require 
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ga 2=% _ 2-9-4) _ 
1+4 P+4 
Using the result of part (i) once more, the steady-state 
solution is 
a = Re((—2— i)e”"*) 
= —2cos 2t + sin 2t. 
To find the amplitude A and phase ¢ of this vibration, we 
express the phasor Z in its exponential form Z = Ae'*, This 
gives 
A= V(-2? +2 = V5 ~224m 
and  @ = —arceos(—2/V5) ~ —2.68 rad, 
so the amplitude and phase are respectively 2.24 m and 
—2.68 rad. 


2h 


5. We seek the phasor Z for the steady-state solution of 


Putting x = Re(Ze'*) leads to + = Re(iwZe™*) and 
# = Re(—w*Ze""), so the phasor Z is given by 

(-w? + awow + wi)Z = P/m 
ene 
wR — w? + awow 

P(wa — w? — i2awow) 

Hence the amplitude of the steady-state solution is 

A=lzi= P\/(e3 — 0)? + da2aRu? 

n((w3 — w*)? + dotweut) 
ps Pim 


and the phase is 


B+ 2ourd +uge =~ coaut = Re (Ze 
= m 


or Z= 


a 
Ved — 0)? + datput 


6. Dividing each coefficient in the differential equation by 
the coefficient of % (so that the coefficient of ¢ becomes 
unity), we obtain 
4+ fa+ 250 
We now compare the coefficients of this equation with those 
of Equation (10) in Section 2, from which it follows that 
2ows = %, we =25, P/m= 8. 
We also deduce from the two right-hand sides that w 
(i) From above, the undamped angular frequency is wo = 5, 
so that the damping ratio is 


@ = Arg(Z) = —arccos ( 


2 
$ cos 6t. 


(ii) The amplitude of the steady-state solution is 
a 
V5 = 0)? + doug? 
‘a 2.25 
V(25 — 36)? + (2.25 x 6)? 
2.25 


= ~ 0.129m, 
121 + 182.25 


and the phase is 


we - w? 
$= =arecos | ————2 
Ved — w)? + dot gu? 
= 
= —arecos ( 
(yates + a) 
~~-2.255rad (or —129.2°). 


‘The angular frequency of the steady-state solution is 
w=G6rads"'. 


A= 


7. From Solution 6 above, we have 
Pe k 


==, ==08=25, A~0.129. 
m 4° om 
Hence k/P = 100/9 and 
Ak’ _ 0,129 x 100 
SS = ie, 


8. The relative error will be less than 10% provided that M 
is between 0.9 and 1.1, which from the graph for a = 0.6 
appears to be the case for the approximate range 

0 < w/wo < 0.9. Since 0.9w) = 0.9 x 20 = 18, the 
corresponding frequency range (for w) is between 0 and 
18rads~'. 


9, We use the values m = 800, k = 5 x 10*, r=9 x 10° 
and P = 800 given in Exercise 2, together with the value 
w ~ 15.7 obtained in Solution 2(i). The undamped angular 
frequency is 


=. | ea, /E ROY Pm 
#0 = 4) =) ago = Tl rads, 


and the damping ratio is 


r 9 x 10% 
6 OT: 
° aVink 2800 x 5 x 10" 


{Note that this value of a is just above 1/2, so that there is 
no input frequency which makes the car body resonate. 

The relevant value of w/a is approximately #27 ~ 2.0. The 
corresponding value of M from Figure 5(a) in Section 2 is 


approximately 0.25. Since M = Ak/P, we then have 
A= MP _ 0.25 x 800 
5x 10° 
‘The estimated value of @ from Figure 5(b) is 
= -0.75m ~ —2.36 rad. 


‘The estimated values of A and @ here are both close to the 
values calculated in Exercise 2. 


=4x 107. 


10. With p = w/w, Equation (13) of Section 2 becomes 
1 


M= 


Differentiation with respect to p gives 
dM _ ~3[2(1—p*)(—2p) + 80*p] 
so that dM/dp = 0 (for points on the graph with horizontal 
slope) when 
— 4p(1 — p*) + 8a*p =0 
or p(p? —14 2a") =0. 
The roots of this cubic equation are 


p=0 and p=+y/1—20?. 
(i) We know that p = w/w» is non-negative. For a < 1/¥2, 
the non-negative roots are p = 0 and p= 1 — 202, and the 
graphs in Figure 5(a) of Section 2 indicate that the second of 
these gives a maximum value for M. (This can be confirmed 
by examining the sign of dM/dp to the left and to the right 
of p= V1— 2a.) 
(ii) For @ = 1/¥2, the root p = 0 is repeated, and for 
a > 1/V2, the expression V1 — 2a? gives complex roots. 
Figure 5(a) shows that in this case M is a decreasing 
function of p for p > 0, so that p = 0 gives the maximum 
value for M. (It can also be shown from above that 
dM/dp <0 for p > 0 and a > 1/¥2.) 
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Solutions to the exercises in Section 3 


1. 


Figure 1 


Tn deriving the equation of motion we assume that the spring 
is extended and that the dashpot’s length is increasing. The 
corresponding force directions are as shown above. 
The extension of the spring is « — lo, and the rate of increase 
of the dashpot’s length is #. The spring and dashpot forces 
therefore have respective magnitudes 

T=k(x-l) and R=rz. 
The only other force acting on the particle is the 
gravitational force, of magnitude mg vertically downwards, 
Newton's second law therefore gives 
mg—-T—-R 
= mg —k(z —lo) — rz 


me 


or mi+rt+kx =mg+klo. 

This equation of motion is identical to Equation (3) in 
Section 1, as expected (since the model considered is the 
same in the two cases), 


ri 
ty 
————— 
NN 
—I*| Post 


We assume, as usual, that the spring is in tension and that 
the dashpot is lengthening. The corresponding force diagram 
is shown above. Note that Pcoswt is the x-component of the 
external force, so that this force has the direction shown only 
when P coswt is positive (the same applies to the diagram 
given in the question), 

The spring force has magnitude T = k(x — Ig) and the 
dashpot force has magnitude R = ré, so that the required 
equation of motion is 

Pcoswt — ri — k(x — lp) 

or mi +ri+ ka = klo + Pcoswt. 


Figure 2 


mt 


3. We seek the phasor Z = Ae“®, where x = Re(Ze'*") is 
the steady-state solution of the differential equation 

m+ ra + kx = Ao(kcoswt — rw sinwt). 
Substituting for x and its derivatives, this equation can be 
written as 

Re|(—mw? + irw + k)Ze""| = Re[Ao(k + irw)e™], 
where the result of Exercise 4(i) of Section 2 has been used 
on the right-hand side. This equation will be satisfied if 


(—mw? + irw + k)Z = Ao(k + irw) 
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or 
z= —Aolk + irw) 
(k= mo*) + ira 
_ Ao(k + irw)[(k — mw?) — irw] 
(R= ma P+ 
Ao ([k(A — mus) + r2w*] — imrs*) 
i (h— mu?) + ra? ; 
This is of the form Z =u —iv, where 
_ Aolke(k — mu) + 120%] 
OS (k= mut)? + tae 


and 
ax Aomrw* 
~ (R= ma??? + rw? 


so that the amplitude is 
A=|2|= fe+u? 
_ Ao y/[e(k = ma?) + 2a? + (marty 


(b= mt? Pa 


v 


and the phase is 
= —arecos(u/|Z|) 
K(k — mu?) + v2? 
= —arccos 5 
[he(he — mw?) + r4w?]? + (mrw)? 


4. (i) When the particle is in equilibrium, the dashpot 

length is constant and so the dashpot exerts no force. The 
force diagram is therefore as shown below, where T is the 

magnitude of the spring force. 


" 


Figure 3 


In equilibrium these two forces must balance, so that 
T = mg. Now if the spring has length 1; when the particle is 
in equilibrium then its extension is /; —/o, giving 
T =k(h —lo). Hence 
k(L, — lo) = mg, 
from which the equilibrium length of the spring is 
hy = lo + mg/k. 
(This is the same argument as that required to answer 
Exercise 1 of Section 1.) 


(ii) 2 


Figure 4 () 
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Figure 4(a) shows the system when displaced from its rest 
position. The displacements z of the particle and y of the 
forcing point are as described in the question. Figure 4(b) 
shows the forces which act on the particle when the spring is 
extended (as it must always be if a rope is being modelled) 
and the particle is moving upwards (that is, the dashpot is 
contracting). The magnitudes of the spring and dashpot 
forces are denoted by T and R respectively. 
The extension of the spring is 

lo+mg/k+y—2-lo, 
so that 

T =k(mg/k+y—-). 
Since # > 0, we also have R = ri. Application of Newton's 
second law gives 

mi=T-—R-—mg 

= k(mg/k + y —2) - rt — mg, 

or mé+ré+ hx =ky. 
This is the required equation of motion. 
(iii) Putting y = Ao coswt in the result of part (ii) gives the 
equation of motion 

mi + ra + kx = kAgcoswt. 
This is the same as Equation (1) of Section 2, except for the 
presence of kA rather than P on the right-hand side, The 
graphs in Figure 5 of Section 2 can therefore be applied to 
this situation where, with kAo in place of P, we have 
M = A/Ao. Note that M here is the ratio of the 
displacement amplitudes of the particle (output) and of the 
upper end of the spring (input). (This quantity is called the 
amplitude ratio of the system.) 
From the values given in the question, we obtain 

wo = Vk/m = V50 = 7.07 rads, 

w=2x 20 > 12.6rads"', 
and hence w/w» ~ 1.78. Corresponding to this value of w/wo, 
the graph for a = 0.2 in Figure 5(a) of Section 2 gives 
M = 0.45, from which 

A= 0A5Ao = 0.45 x 0.05 ~ 0,023 


‘The amplitude of the particle's steady-state oscillations is 
therefore about 0,023 m, or 2.3 cm. 


[We may check that the spring representing the rope is in 
tension throughout the motion by showing that 7’ (as defined 
above) is always positive. Since |x| < A and |y| < Ao, we 
have 

T = mg + k(y — 2) > mg — k(Ao + A). 
The value of the right-hand expression is approximately 

4 x 9.8 — 200(0.05 + 0.023) ~ 24.6N, 


so the spring is always in tension.) 
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Solutions to the exercises in Section 4 


‘The model predicts that a graph of 

log, {(Ao + Ad)/(An + A’,)] against n should be a straight 
line through the origin (with slope pr). The given data have 
heen plotted on such a graph above. There is no straight line 
through the origin which provides a convincing fit to these 
data points, so this test indicates that linear damping is not a 
satisfactory modelling assumption for a water-filled dashpot, 


2. (i) Assume that the spring is extended and that the 

dashpot is lengthening. The spring force is then of magnitude 

T = k(x — y — lo) and the dashpot force is of magnitude 

R= ra, both being directed downwards. The only other 

force acting is that due to gravity, which is of magnitude mg 

directed downwards. Newton's second law therefore gives 
mé = —-T —R—mg 


= =k(x -—y — lo) — rt — mg, 
or mié+ra+ kx = klo — mg + ky, 
which is Equation (1) of the text. 
(ii) Equation (2) of the text is 
mi + ré + ka = k(lo + yo — mg/k) + kAp coswt. 
The steady-state solution is 
© = lo + yo —mg/k + Acos(wt + 4), 
where the sinusoidal term is the steady-state solution of the 
equation 
ma + ra + kx = kA coswt. 
We now proceed much as in Solution 4(iii) of Section 3. The 
last equation is the same as Equation (1) of Section 2, with P 
replaced by kAg. The analysis of Subsection 2.3 can therefore 
be applied, leading to Equation (13) of Section 2, which was 
————— 
VU = (w/wo))? + 402 (w/w)? 
Here M is equal to A/Ao, since P has been replaced by kA. 
This is the expression for the amplitude ratio given in 
Equation (3) of Section 4. 


(iii) From the result of Exercise 10 of Section 2, there will be 
no resonance provided that the damping ratio a is greater 
than or equal to 1/72. 
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Figure 2 


We assume that the spring is extended and that the dashpot 
is lengthening. Then the spring and dashpot forces are as 
indicated in the diagram. The lengths of the spring and 
dashpot are both x — y. Hence the extension of the spring is 
a2 —y—lo, and the rate of increase of the dashpot’s length is 
&—j. The magnitudes of the spring and dashpot forces are 
therefore 


T=ka-y-lo) and R=r(z-y) 
respectively. So the equation of motion of the particle is 

mé = -T —R-mg 

= —k(a — y—lo) — r(é- 9) — mg 

or mi+re+ kx =klo —mg+ryt ky. 
(ii) Since = y, we have 

i=a-y9 and F=7-j. 
‘The equation of motion becomes 

m(i+ fj) +r(2+y) +k(z+y) =klo — mg + ryt ky 
or mi+ri+kz =klo — mg — mj. 


4. The equation of motion of the tone-arm in terms of 
2 =x~—y is given by the last equation in Solution 3 above. 
With y = Ao coswt we have jj = —w* Ap coswt, giving 
mi + ri + kz = klo — mg + mw* Ap coswt. 
The steady-state solution of this equation is 
= = Ip — mg/k + Acos(wt + 4), 
where the sinusoidal term is the steady-state solution of 
mE +ri+kz = mw" Ao coswt, a) 
This is, once again, the same as Equation (1) of Section 2, 
but with mw*Ap in place of P. As before we may therefore 
conclude that Equation (13) of Section 2 holds, which is 
1 
V0= @/o0)?)? + 403(w/o)? 
Here M = Ak/P (from Section 2) is replaced by 
Ak _ uA 
muAg  W?Ag* 


since wo = /k/m. Hence we have the amplitude ratio 


Ae (w/wo)? 
Ao J — w/uo)?)? + 4a7(w/an)? 
as required. 


As an alternative, we show below how the phasor method 
can be used to derive this expression directly, without 
reference to Section 2. Starting from Equation (1) above, we 
divide through by m and then re-express the result in terms 
of the undamped angular frequency wo = 
damping ratio a = r/(2Vmk). This gives 
= + 2aweé + waz = w* Ao cost. 
We seek a solution to this equation of the form 
2 = Re(Ze'“"), where Z is a complex constant. Since 
2 = Re(iwZe“*) and = = Re(—w*Ze""), substitution into 
the differential equation produces 
Re[(—w? + Qiawow + ws )Ze™*] = Re(w? Age™*). 
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This is satisfied provided that 
(-w? + 2iawow + 03)Z = w* Ao, 
for which we require 
wo 
w2 — w? + Ziowow 
= (w/o)? Ao 
~ T= (e/a)? + 2iaw uy" 
The amplitude of these steady-state vibrations is 
(w/wo)? Ao 
VG = (w/ic0)?)? + 403(0/0)? 


giving the required amplitude ratio as before. 


Z= 


A=(Z|= 


Solutions to the exercises in Section 5 


1. (i) We model the object by a particle of mass m, and 
assume that the spring and dashpots are perfect, with the 
spring having stiffness k and each dashpot having dashpot 
constant r. The particle displacement x is measured 
downwards from the equilibrium position. This model 
system is shown on the left below. 


a 


Figure 1 


Assuming that the spring is extended and the dashpots are 
contracting, the forces acting on the particle are as shown on 
the right of the diagram above. The magnitudes of the 
spring and dashpot forces are respectively 

T=mg+ker 
(since T = mg when the particle is in equilibrium) and 

R=rlé|=re. 
Newton's second law then gives 

mé=mg-T-R-R 

= mg — (mg + kx) — 2ré 

or mé + 2ré + ka =0. 
This can be written as 

#+2awot + wor = 0, 
where wy = \/k/m is the undamped angular frequency and 
a=r/Vmk is the total effective damping ratio (that is, the 
ratio obtained by regarding all the damping as being 
produced by a single dashpot with dashpot constant 2r). In 
the equilibrium position, a force of magnitude mg produces a 
spring extension of 19.6cm = 0.196 m, so that 


mg = 0.196k 
= 50. 


i g 
oF HO = Tn ~ 0.196 
Also a = 1 is given, so the equation of motion is 
#+2V50z + 502 = 0. 
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(ii) Removing one dashpot has the effect of halving the 
damping. The new equation of motion is therefore 


+ V50E + 50x = 0, 
where the damping ratio is now a = r/(2Vmk) = 0.5. From 


Equations (7) and (8) of Section 1, the general solution of 
this equation is 


x= Ae~* cos(Qt + ¢). (1) 
where 
p=, n= Ynk=F 
ae 
and A, @ are arbitrary constants. From w) = 50 and a = 0.5, 
we have 
k=50m and r=Vmk=V50m, 
giving 
p=VI25 and 2= 4200-50 = V375. (2) 
The period of damped vibrations is 
Qn an 
7 = y= Japp = 18 


If the particle starts from rest with the spring at its natural 
length then, since the origin is at the particle’s equilibrium 
position, the initial conditions are 


x(0)=—0.196 and #(0)=0, 


Putting ¢ = 0 into Equation (1) and applying the first of 
these conditions produces 


Acos ¢ = —0.196. 
‘The maximum subsequent spring extension occurs after half 
a period at time t = 37 = /2, when (t) = 0 once more. 
The corresponding value of x is (using Equations (2)) 
2(x/2) = Ae~?"!" cos(x + 6) 
= —Ae?"!" cos 
= 0.196e~"/¥5 ~ 0,032. 
The spring extension is then approximately 
0.196 + 0.032 = 0.228 m, 
or approximately 22.8 cm. 


2. The equation of motion for Exercise 1(ii) was 

mi +ré + ke =0, 
and with the additional downward magnetic force, this 
becomes 

mi +ra + ka = Py + Poosut. 
The steady-state solution is 

2 = Po/k + Acos(wt + ¢), 
where the sinusoidal term is the steady-state solution of the 
equation 


mi +ri+ kx = Pcoswt, 
that is, 


#+2aun# +ub2 = © cosut, (3) 


where wo = \/k/m and a = r/(2Vmk). Equation (3) is 
identical to Equation (10) of Section 2, and the 
corresponding amplitude of forced vibrations is given by 
Equation (11) of Section 2 as 


Pim 
J(u? = 0)? + datwiu? 


Now a = 0.5 from Exercise 1(ii), so the amplitude is 


we P/(rnsi) 
VG= feo)" + 70)? 
d 


| Va-r P+ 
where p = w/w and d = P/(mw%) = P/k is the distance by 
which a downward force of constant magnitude P extends 
the spring. 
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3. (i) Let x be the upward displacement of the particle, 
measured from the fixed point about which the sinusoidal 
motion y = Ag cos wt of the base of the box takes place. 
Then we have «= y + = +d, where d is the length of the 
dashpot in the equilibrium position. Assuming that the 
spring is extended and the dashpot lengthening, the forces 
acting on the particle are as shown in the diagram below. 


ae ee 


Figure 2 


‘The spring and dashpot forces are respectively 

T=mg—kz 
(since T = mg when the particle is in equilibrium) and 

R 2 
Newton's second law therefore gives 

méi=T-—mg-—R 

= (mg — kz) — mg — ri. 

On putting x = y + = +d we have # = j + 2, so the equation 
of motion relative to the box is 


mi+ri+k mip 
or mi+ré+kz = mw" Ap coswt. 


(ii) The last equation is a copy of Equation (1) of Section 2, 
with mw® Ag in place of P, so by Equations (11) and (12) of 
that section we have the steady-state solution 
== Acos(wt +), 
where 
w* Ao 


(a — wy? + detusut 


ee ~ ae ( 


we —w? 
Ved — 0? + doPadu? 


MST204 8 Solutions 


(iii) If wo /w and a are both small, then 
= Aa ee 
= = Ao 
V ((wo/w)? = 1)? + 40?(wo /w)® 


and 


see (joreeee) 
V (w/w)? — 1)? + 4a (wo/w)? 
~ —arccos(—1) = —7. 
Hence we have 
2 = Acos(wt + o) 
~ Ap cos(wt — 7) = —Ag coswt 
=y. 


or 


4. (i) If w/w» and a are both small, then 
(w/wo)?Ao ~ wt Ao 


Vi = who)? +407 (woo 


and 


| d= (w/er)? ) 


V (1 = (w/w0)?)? + 403 (w/w)? 
= —arecos(1) = 0. 
Hence we have 


Acos(wt + ) 
~ (w/wo)* Ao cos wt 


or 2% —G§/ug. 
Thus = is proportional to jj, and the constant of 
proportionality —1/a@ is independent of w. 


(ii) The approximate value of acceleration as indicated by 
the instrument is 


wiz = —weAcos(wt + 4), 


whose amplitude is wj A. The actual value of the acceleration 
is 


i = —w* Ao coswt, 


with amplitude w* Ap. Using the expression above for A and 
putting a = w/wo = 0.1, the ratio of these amplitudes is 


weA 7 1 
WA (1 = (w/wo)?)? + 40 (w/w)? 
i 


~ Jaap? +4x 018 


= 1.01. 


